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Abstract 

Let $R be a complete discrete valuation ring, S = £H[[u]] and n a positive integer. The 
aim of this paper is to explain how to compute efficiently usual operations such as sum and 
intersection of sub-S-modules of S d . As S is not principal, it is not possible to have a uniform 
bound on the number of generators of the modules resulting of these operations. We explain 
how to mitigate this problem, following an idea of Iwasawa, by computing an approximation 
of the result of these operations up to a quasi-isomorphism. In the course of the analysis 
of the p-adic and w-adic precisions of the computations, we have to introduce more general 
coefficient rings that may be interesting for their own sake. Being able to perform linear 
algebra operations modulo quasi-isomorphism with S-modules has applications in Iwasawa 
theory and p-adic Hodge theory. 
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1 Introduction 

Let 91 be a complete discrete valuation ring (see ij2.ll for a reminder of the definition) whose 
valuation is denoted by . Let K denote its fraction field with valuation Vk and 7r be a uniformizer 
of $H. We set S = 9t[[u]]; it is the ring of formal series over *H. Our aim is to provide efficient 
algorithms to deal with finitely generated modules over S. Since, we can always represent a torsion 
module as the quotient of two torsion-free modules, we shall focus on torsion-free modules. 

Any finitely generated torsion-free 5-module j$ can be considered as a submodule of S d for d 
big enough. As a consequence, we can represent ^# by a matrix whose columns are the coefficients 
of generators of ^ in the canonical basis of S d . Thus we can reformulate our problem as follows: 
given Mi and two matrices representing respectively the S'-modules J£\ and embedded 
in S d , give algorithms to compute a matrix representing n ^#2 or + ./#2- We would like 
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also to be able to check membership, equality of sub-S-modules, inclusions, etc. As S is not a 
principal ideal domain, in order to control the number of generators of the sub-S'-modules of S d , 
we propose, following an idea of Iwasawa, to compute approximations of the submodules resulting 
of aforementioned operations in the following sense: we say that a morphism ^#1 — ¥ is a quasi- 
isomorphism if its kernel and co-kernel have both finite length, and we want to make computations 
modulo quasi-isomorphisms. We propose two different approaches, each of them having its own 
advantages and disadvantages. 

First, we notice that there exists a correspondence between the set of classes of modules modulo 
quasi-isomorphism and modules over the rings S n and S u defined respectively as the localization 
of S with respect to tt and the completion of the localization of S with respect to u. For A = 
S^, S u , let Free^ be the set of free sub-A-modules of A d , and denote by Mods/ qis the set of quasi- 
isomorphism classes of sub-S'-modules of S d , there is an injective morphism ty' : Modg/ qis — > 
Free^ x Free^, Jl h-> (Jt ®s S ff , J% ®s S u ), where ^# is any representative in the class The 
image of <£' can be precisely characterized (see Theorem 11.11 below 1 ). Using this correspondence, 
operations with modules with coefficients in S reduces to the computation with modules over 
and S u . As these two last rings are Euclidean, there exist classical canonical representations and 
algorithms to manipulate modules over these rings. 

A second approach consists in finding a canonical representative in a class of modules modulo 
quasi-isomophism which is amenable to computations. Such a representative is provided by the 
maximal module of a S-modulc ./#. It can be defined as the unique free module in the class of 
quasi-isomorphism of ^ . We present an algorithm to compute the maximal module associated 
to a sub- S- module of S d which is inspired by a construction of Cohen, presented in (TU1 p. 131], 
to obtain a classification up to quasi-isomorphism of finitely generated S-modulcs. We can then 
compose this algorithm with algorithms to compute basic operations on free modules in order to 
compute with representatives up to quasi-isomorphisms. 

In order to obtain real algorithms (i.e. something computable by a Turing machine) we have 
to consider the fact that elements of S and its localized are not finite. In this paper we consider 
an approach in two steps in order to solve this problem. First, we give the ability to Turing 
machines, to manipulate, by the way of oracles, elements of S, S^, S u . More precisely, we suppose 
given oracles able to store elements of the base ring, compute valuation, multiplication, addition, 
inversion, and Euclidean division. We express the complexity of an algorithm with oracle by the 
number of calls to the oracles to compute ring operations. Once we have well defined algorithm 
with oracles to compute with modules, we study in a second time the problem of turning them 
into real algorithms. 

Much in the same way as for floating point arithmetic, the actual computations with modules 
with coefficients in S are done with approximations up to certain 7r-adic and it-adic precisions. It 
is important to ensure that the (truncated) outputs of our algorithms are correct which means 
that they do not depend on the tt or u powers of the input that we have forgotten. In order 
to deal with this precision analysis, it is convenient to consider a generalisation of the family 
of ring coefficients S. Namely, given a, f3 relatively prime integers, we write v = j3/a and set 
Si, = {J2 a i ul G ^[[w]]|vR-(ai) + vi > 0, Vi G N}. We have So = S. In this paper, we develop 
a theory of S^-modules which encompass modules over S and use it in order to obtain algorithm 
with complexity bounds and proof of correctness. 

More precisely, we generalize the definition of a maximal module for finitely generated torsion- 
free Sj,-modules. Denote by Max^ the set of maximal sub-S„-modules of S d . We prove the 
following theorem (see Theorem 13. 12)) . which generalize the above mentioned decomposition: 

Theorem 1.1. The natural map 



is injective and its image consists of pairs (A, B) such that A and B generate the same S ' -vector 
space in S d . If a pair (A, B) satisfies this condition, its unique preimage under 'F is given by Ap\B. 

In the theorem, S is a field containing S v and its localized S v ^ and S„, M which is precisely defined 
in Section [2~2l We give an algorithm with oracles to compute the maximal module associated to 
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a finitely generated torsion-free S^-module. In general, it is not true that the maximal module 
of a torsion-free S^-module is free, although this property holds when v = 0. Nonetheless, by 
using the theory of continued fraction, it is possible to obtain a tight upper bound on the number 
of generators of a maximal module embedded in S%. If v is rational, it admits a unique finite 
development as a continued fraction that we denote by [ao; <Xi, . . . , a n ] (here, we suppose that 
a n ^ 1). We can prove the following (see Theorem 13. 32)) : 

Theorem 1.2. Let v = [ao;a±, . . . ,a n ]. Let ^ be a sub-S^ -module of S^. Then Max(^#) is 
generated by at most d • (2 + X)I={ 2 ^ a 2i) elements. 

We provide some simple examples to show that a lot of basic operations that we need in order 
to compute with modules over S v , such as the computation of the Gauss valuation, are not stable. 
This means that, in general, the computation with approximations of the input data does not yield 
approximation of the result. This is where it becomes interesting to use the possibility to change 
the slope v of the base ring S v . In the context of our computation, a bigger slope plays the role of 
a loss of precision in the computation of an approximation of a module over S v . In this direction, 
we can prove the following theorem (see Theorem 14.71 for a precise statement): 

Theorem 1.3. Let jM\ and ^#2 be two finitely generated sub-S^ -modules of Sf, such that ^2 C 

1 / 7t c ^#l for a positive integer c. Let Mi and Mi be the matrices with coefficients in S v of generators 
of and -#2 in the canonical basis of S„. Suppose we are given approximations M[ and M^ of 
M\ and M2 respectively. Then, for a well chosen v' > v, there exists a polynomial time algorithm in 
the length of the representation of M[ and M£ to compute a matrix M% which is an approximation 
of the maximal module associated to (y^i C8s„ S v >) + (.-#2 <S>s„ S v >). 

The organisation of the paper is as follows: in fj2j we introduce the rings S Ul and their basic 
arithmetic and analytic properties. In SJ21 we generalize some classical results of Iwasawa to the 
case of finitely generated ^-modules and then give an algorithm with oracle to compute the 
maximal module associated to a torsion-free ^-module and obtain an upper bound on the number 
of generators of a maximal module. Note that iJUand <J21 we only describe algorithms with oracles. 
In 21 we study the problem of p-adic and u-adic precisions and turn the algorithms with oracles 
obtained in the previous sections into real algorithms. 

2 Arithmetic of the rings S u 

In order to compute with modules over S v we first have to study the basic arithmetic properties 
of their base ring. In this section, we show that its localized with respect to u a /ir^ and tt becomes 
Euclidean. We provide algorithms with oracles to compute the Euclidean division in these rings 
which will be very useful for our purpose along with their complexity expressed in term of the 
number of ring operations. They will be turned into real algorithms (i.e. working on a real Turing 
machine) in ij4] where we study the problem of precision of computation in the rings S v . 

2.1 Notations 

We fix the notations for the rest of the paper. Let £H be a ring equipped with a discrete valuation 
that is a map w<k : SH —¥ N>o U {+00} satisfying the following conditions: 

• for all x £ v^\(x) = +00 if and only if x = 0; 

• for all x G £H, v<r(x) — if and only if x is invertiblc; 

• for all x.y <G $H, v^(xy) = v<x(x) +v SK (y); 

• for all x, y G v<ji(x + y) > ia.m.(v<x_(x) , u<r(j/)). 

Let a be a fixed real number in (0, 1). One can define a distance d on *H by the formula d(x, y) = 
a vm(x~y) [x^y g 5H) where we use the convention that a + °° — 0. For the rest of the paper, we 
assume that D\ is complete with respect to d. We recall that D\ is a local ring whose maximal ideal 
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is DJl = {x G > 0}. Up to renormalizing v<r, it is safe to assume that it is surjective, what 

we do. We denote by it a uniformizer of Dt, that is an element of 1H whose valuation is 1. Every 
element x in d\ can then be written x = tt t u where r = vsr(x) and u G 9i is invertible. Here are 
several classical examples of such rings D\: 

• the ring 1 V of p-adic integers equipped with the usual p-adic valuation; 

• more generally, the ring of integers of any finite extension of Q p ; 

• for any field k, the ring k[[u]] of formal power series with coefficients in k. 

We now go back to a general £R. It follows easily from the definition that the field of fractions of 
£H is just Dt[l/7r]. Let's denote it by K and set S — 9t[[u]], the ring of formal series over 9i. The 
valuation v?x can be extended uniquely to a valuation vk on K. 

2.2 Definition and first properties of S v 

Denote by if [[it]] the power series ring with coefficients in if. It is classical to define the Gauss 
valuation of an element ^aiU 1 G K[[u]] as the smallest vk((h) if it exists. The ring of elements 
of if [[u]] with non negative Gauss valuation is nothing but In this section, we are going to 

consider more generally a family of valuations parametrized by a slope v G Q so as to define the 
subring of if [[u]] of elements with positive valuation. 

Definition 2.1. Let i/ € Q. We define the Gauss valuation v v : K[[u]] — s- Q U {+oo, — oo} by 
v v (x) = +oo if x = 0, a i ul ) = mm { u -R"( a i) + vi, i £ N} if x ^ and this minimum exists and 

v v (x) — — oo otherwise. The Weierstrass degree of x denoted deg^(x) is given by deg^(0) = — oo, 
deg^(x) = mm{i\vK(ai) + vi = v v (x)\ if v v (x) ^ — oo and deg'{ v (x) — +oo otherwise. When no 
confusion is possible, we will use the notation deg w instead o/deg^/. 
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Figure 1: The Gauss valuation of 7r 2 • u 4 with v = 1/3 is 10/3. 

The following lemma gives some basic properties of v v and deg w . In particular, it shows that 
v v has the usual properties of a valuation: 

Lemma 2.2. For x,y G if[[u]] we have: 

1. v v {x) = +oo if and only if x = 0; 

2. v v (x ■ y) = v v {x) + v u (y); 

3. v v (x + y) >mm{v v (x),v v {y)). 

Moreover for all x, y G if [[u]] with finite Gauss valuation, deg w (x.y) — deg w (x) + deg w (y). 

Proof. To prove 2., we first suppose that x = ^a^u 1 and y — ^btu 1 have finite valuation. Let 
z = x ■ y = ^CiU 1 . We have vn(ci) + vi = VK(Ylj=o a i ' bi-j) + vi > minj{«x(oj) + v ■ j + 
vic(bi-j) + v ' [i > v v{x) +v v {y). Moreover, by taking i = deg w (x) + deg w (y) in the previous 
computation, we obtain that v k {cd C g w (x)+dcg w (y)) — v v (x) + v lJ (y). If v v (x) — — oo and y ^ 0, we 
can apply the previous result to the series obtained by truncating a; up to a certain power to show 
that v v (x ■ y) — — oo. The proof of the rest of the lemma is left to the reader. □ 
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We let S v = {x £ if [[tt]]|?;„(x) > 0}. By definition, an element x £ S u can we written as a 
series 

X = dill 1 , 

ieN 

where a, £ if and vxiai) > —vi. 

Remark 2.3. i£ is c/ear i/iai S v is complete for the valuation v v . Nonetheless, the ring S u is not 
a valuation ring. In fact, although v u (u a /-k^) — for v ^ (resp. v v {u) = for v = 0), u a /ir 13 
(resp. u) is not invertible in S v . 

We let 

S u ,tt = S v [l/ir] = { aiU l ,ai £ if such thai Vx(ai) + Abounded below}. 

ieN 

In the same way, it is clear that one can extend the v v valuation of S v over S v \n& /u a ] and we 

let Sjs. u = Sj/ [k^ /u a ] where the hat stands for the completion of S^tt 13 /u a ] with respect to the 
topology defined by v v . 
Put in another way, 

S V u = { / o-iU 1 , at £ S v and lim Vfc(ai) + vi = +00} . 

We moreover define 

S ={ Voiii^fli G if Vk{o,z) + vi bounded below and lim Vxio-i) + vi = +00}. 

^ — ' i^ — oo 

We have the following commutative diagram of inclusions: 



(1) 



As S v>7r is a subring of if [[«]], it is equipped with the v v valuation and the Weierstrass degree 
associated to v v . Moreover, one can extend, in an obvious manner, the definition of v v and the 
Weierstrass degree for S VjU and S. 

We can interpret the ring S v in terms of the analytic functions on the 7r-adic disc. In order 
to explain this, for v £ Q, we consider the open disk D u = {x £ K\vk(x) > v}. Denote by & v 
the ring of convergent series G v = {X)ieN a i u% \ a i ^ if, hniinfi^ +00 VK ^ a ^ > —i/} in the disk D u . 
It is clear that S v ^ is exactly the set {/ £ if [[u]] \ Vk {f{x)) bounded below on D v } and S v can be 
described as {/ £ if [[u}) \ vk{J{x)) bounded below by onD„}. Thus, there are obvious inclusions 
S v C S„ )7r C & v but one should beware of the fact that the last inclusion is strict. Indeed for 
instance, for £H = Z p , v = the series X)i>o T wri ich defines the function log(l — u) is convergent 
in the unity disk but is obviously not in So l7r since v^(l/i) has no lower bound. Assuming that v is 
rational (what we do), the following proposition gives another characterisation of elements of G v 
that lies in S„ )W . 

Proposition 2.4. An element x £ G v is in S v ^ if and only if x has only a finite number of zeros 
in O v . 

Proof. Let x £ G v . The number of zeros of x £ D v is equal to the length of the interval above 
which the Newton polygon of x has a slope < —v. If this length is finite, it is clear that v p (di) is 
bounded below by a line of the form — vi + c with c a constant and as a consequence is an element 
of S v ^. 

Conversely, suppose that x £ S v ^. This means that v p (ai) + vi is bounded below and is 
contained in Z+vZ which is a discrete subgroup of R (as v is rational). Thus, the set {v p (ai)+vi, i £ 
TV} reaches a minimum for a certain index iq. This means that for all i > io, the slope of the Newton 
polygon of x is greater than — v and x has a finite number of zeros in D v . □ 
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We end up this section, by remarking that up to an extension of the base ring U\ all the S^s 
are isomorphic to a Sq. Indeed, write v = fj/a with a, (3 relatively prime numbers and let w, in an 
algebraic closure of K, be such that vo a = ir. Let 9V = £H|W], K' be the fraction field of 5H' (and a 
finite extension of K). The valuation on £R extends uniquely on 9V by setting VK'ipz) — I /a. For 
/-t = 0, v, let = S 1 ^ 0<h 9^' ■ The valuation defines a Gauss valuation on S^' that we denote 
also by v^. 

Lemma 2.5. Keeping the notations from above, the morphism of ring p : So —> S' v , defined by 
p(l) — 1 and p(u) = -T- is an isomorphism. Moreover, if x 6 So' we /icrae WoOO = ^(^(a;)) and 
deg^(x) = de$p(p(x)). 

Proof. By definition, S' v — {^2 aiU i \vK' (a>i) + vi > 0} = {^2,ai{u/w^y\vK'{<ii) > 0} from which it 
is clear that p is an isomorphism. The rest of the lemma is an easy verification. □ 



2.3 Division in S u 

The Weierstrass degree allows us to describe an Euclidean division in S u . Although, the existence 
of such a division is classical (see for instance [TU]) at least over So = £K[[u]], we give here a proof 
for all v which provides an algorithm with oracles to compute the Euclidean division. 
In order to study divisibility in S„, we have a first result: 

Lemma 2.6. Let x.z £ S v . We suppose that deg w {x) = then there exists y £ S v such that 
x.y = z if and only if v v (x) < v v (z). 

Proof. We suppose that deg w (x) = 0. If there exists y £ S v such that x ■ y = z then clearly 
v u (z) < v u (z). Reciprocally, we suppose that v v (x) < v u (z). Write x = J2i£fi a i u '' an( I z = 
SieN c i u% ■ Since oq is invertible in K there exists y £ ^[[u]] such that x.y = z. We have to 
prove that v v (y) > 0. For this, write y = Y^ie^hW 1 - We have Wrt(&o) = Vk(cq) — vk(o,o) > 
by hypothesis. Then, for j > 1, we prove by induction that vicipj) + vj > 0. We have bj = 
Oq 1 ■ cj — Oq 1 Yjl=i a i ' bj—i- But vk{olq 1 ■ cj) + uj > v v (z) — v„(x) > because deg w (x) = 0. 
Moreover, for i = 1 .. ,j, vk{o-q 1 • • + KJ = ^A'(a;) + — + VK{bj-i) + v{j — «)• But 

by definition vxifli) + vi ~ v v (x) > and by the induction hypothesis vn(bj-i) + v(j — %) > 0. 
Therefore, vx{bj) + vj > and we are done. □ 

Applying Lemma 12.61 to z = 1, we get 

Corollary 2.7. Let x — a i xl e &v> then x is invertible in S v if and only if deg w {x) = and 

v u (x) = 0. 

We note that the corollary implies that S v is a local ring. Next, we introduce the following 
notations: for x — X^eN a * M * e ana ^ ^ a positive integer, we let Hi(x, d) — X^i>d a i wI an d 
Lo(a;, d) = J2t=o a i ul - It is clear that x — Lo(x, d) + Hi(x, d). 

Proposition 2.8. Let x,y £ S u . Suppose that v v {y) > v u (x) then there exist a unique couple 
(q, r) £ S v x (K[u) fl S„) such that deg(r) < deg w (x) and y — q ■ x + r. 

Proof. First, we prove the existance of (q, r). Let d = deg w (x), we consider the sequences (qt) and 
(ji) defined by qo = and ro — y and 

m(n,d) Hi(r i; d) 
= * + HiM' r ' +1 = n - ffiM ' X - (2) 

We are going to prove by induction that qi and are convergent sequences (for the valuation) 
of elements of S v . Let e = v v {Lo(x, d)) — v v (Si[x, d)) > 0. Our induction hypothesis is that qi and 
Ti are elements of S v , that u w (Hi(r,-, d)) > e • i + f^(Hi(y, d)) and that y — qi ■ x + ri. It is clearly 
true for i = 0. 

By the induction hypothesis, we have v l/ (Ri(r i ,d)) > v v (Hi(y, d)) and by hypothesis v v (Hi(y, d)) > 
Vv{y) > v v {x) — v u (Hi(x, d)) so that v v {Hi(ri, d)) > v u (Hi(x, d)). Applying Lemma \'2.G\ we obtain 
m ( a 'd) e ^ ana * then rj+i £ 5^. Next writing x — Hi(x, d) + Lo(x, d), we get 

r i+ i = Lo(n,d) - ^TT^jr -Lo(x,d). (3) 
Hi(x, d) 



G 



Applying Lemma f2.21 we obtain that u„(Hi(rj+i, d)) > u !/ (Hi(rj, d)) + v v (Lo(x, d)) — v v (Hi(x,d)). 
Using the induction hypothesis, we get that v„(Hi(r,+i, d)) > e ■ (i + 1) + v v (Hi(y,d)), Finally, 
using the hypothesis that y = qi ■ x + ri, we immediately check using ([2]) that y = qi+\ ■ x + r,-+i. 

From the induction, we deduce that qi and r, are convergent sequences of S v for the v v valuation. 
In fact, we have q i+1 - q t = "'^^ so that v u (q i+ i - q,) = v u (Ui(n,d)) - v u (Hi(x,d)) > e ■ i + 
v u (Hi(y, d)) — v u (Hi(x, d)) > e ■ i. The same argument works for r^. Denote by q and r the limits. 
As for all i G N, y — • x + n, we have y = q ■ x + r. Moreover, since Hi(ri, <i) > e • i, we have 
Hi(r, d) — 0, so that r S /sT[m] and deg(r) < deg VK (a;). 

We prove the unicity of (g, r). Let (q 1 , r') G 5*^ x (7\T [u] n 5^) such that y = q' ■ x + r' . Then 
(g — g') • x = r' — r. We have deg w ((q — q') ■ x) = deg w (r' — r) < deg w (x) which is only possible 
if q = q' and r = r' . □ 

From the proof of Proposition 12.81 we deduce Algorithm [1] to compute from the knowledge of 
x, y, the elements q',r' G S u such that v v {q — q') > prec and v v (r — r') > prec. Furthermore, by 
the proof of the proposition, the number of iterations of the while loop is bounded by \prec/e\. 
We deduce that Algorithm [T] needs one inversion and 3 • \prec/e] multiplications in S v . 



Algorithm 1: EuclieanDivision 



input : x,y G S u with v v {y) > v„(x), prec G N 

output :g,reS„ such that y = q ■ x + r and tv(Hi(r, deg w (x))) > prec 

1 q <- 0; 

2 r <- y; 

3 d <r- deg w (x); 

4 while v„(Ri(r, d)) < prec do 

5 



Hi(r,d) 

_ m(r,d) 
r *~ T Hi(x,d) ' 



6 

7 return q, r; 



Now, let x G iSy, following [lOj we say that x is distinguished if tv(x) = 0. With this definition, 
we can state the classical Weierstrass preparation theorem: 

Corollary 2.9 (Weierstrass preparation). Let x G S u be a distinguished element and let d = 
deg w (x). Then we can write x = q ■ h, where q G S v is an invertible element and h G K[u] Pi S„ is 
of the form h — + Y^i=o ^ iU% w ^ n v K(bi) + vi > 0. 

Proof. We first notice that dv is a nonnegative integer. Indeed, it is clearly nonnegative, and 
writing x = J2 a dU d , we have v<j\(o>d) + dv = (since x is assumed to be distinguished) and, 
consequently, dv = —vy^{a c i) G Z. 

By proposition 12.81 there exist q& S v and r G K[u] f] S v such that degr < d and 

u d 

—— — q ■ x + r. 

7T 

Using Lemma \2~2~[ we obtain Vj,(<z) = and deg w (q) = 0. Then, Corollary 12.71 implies that q is 
invertible. To finish the proof it suffices to remark that deg w (-^-rr — r) = d and the result follows 
from the definition of deg w . □ 

Remark 2.10. The previous proposition is closely related to the Proposition \2.4\ since it says that 
an element of & v is in if and only if it can be written as product of a polynomial times a 
function which does not have any zero in D v . 

The following proposition states that the rings S^.^ and S v . u are Euclidean rings and provides 
algorithms with oracles to compute the division. 

Proposition 2.11. The ring S v ^ is Euclidean, the ring S VjU is a discrete valuation ring for the 
valuation v v (and as a consequence is also Euclidean). Moreover, § is afield. 
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Proof. Let x, y £ S„ tW . There exist s.ieff such that ir s x,Tr t y S 5jy and iv(tt* ■ y) > w„(-7r s • sc). 
Applying Proposition 12. 8[ yields q e S v and r e [[u]] H 5^ such that deg(r) < deg w (x) and 
y = 7r s_ * • q ■ x + 7r~* • r and we are done. 

In order to prove that S VyU is a discrete valuation ring, we have to show that the set of invertible 
elements of S v , u is the set of elements x 6 S u , u such that v v (x) — 0. Write v = f3/a, with at,j3 
relatively prime numbers. Let m be the ideal defined by {x 6 S UtU ,v u (x) > 0}, it is clear that 
S V}U /™ is isomorphic to the field k((u a )). As S UyU is complete for the v v valuation, the Hensel lift 
algorithm gives an algorithm with oracles to compute the inverse of an element whose valuation 
is zero. The Algorithm [2] uses a fast Newton iteration to perform this computation modulo m™ at 
the expense of 0(log(n)) multiplications in S u , u . 

Let a; be a non zero element of #, by dividing it by a power of it we can suppose that v u (x) = 
and by using the algorithm with oracle Algorithm we can invert it. □ 



Algorithm 2: Inverse 



input : x € S v ,u such that v u (x) =0,n£N 
output : y € S v , u such that x ■ y = 1 mod m r 

1 if ?i = 1 then 

2 | y -ir- 1/x mod m; 

3 else 

4 y <— Inverse(x, [n/2]); 

5 y <— y + y(l — ay) mod m"; 



Remark 2.12. One can use the usual Euclidean algorithm to compute the Bezout coefficients of 
!,!/£ Sy^. This algorithm outputs g,k,l,m,n e S^.^ such that g is the greatest common divisor 
of x and y, k ■ x + 1 ■ y — g, m ■ x + n ■ y — and k ■ n— I ■ m = 1. It proceeds by using the fact that 
gcd{x, y) — gcd{y, r) where r is the rest of the division of x by y and uses 0(deg w (y)) calls to the 
Euclidean division Algorithm [7J We remark, as the rest of the division of two elements of S v is 
an element of K[u], that starting from the second iteration of this algorithm all the divisions to be 
computed are the usual division between elements of K[u]. Unfortunately, we will see that in fj^J that 
the Euclidean algorithm in general is not stable, so that we might need extra informations, about x 
and y in order to compute an approximation of their gcd from the knowledge of an approximation 
of x and an approximation of y. 



3 Modules over S v 

Let d be a positive integer and fixvgQ. We want to compute with finitely generated torsion free 
S'jy-modules. Any such module ^# can be embedded in S£ for d 6 N and can be represented by 
a matrix with coefficients in S v whose column vectors are the coordinates of generators of ^ in 
the canonical basis of 5^. Indeed, we can always embed ^# is ^# <S>s u Frac(S , !y ) and select a basis 
(ei, . . . , ed) of Jt ®s v Frac(S„) together with an element D <E S v such that the image of j$ in 
®s v Frac(S' t/ ) is contained in the free S^-module generated by the -^.ei's. 

A first problem arises here: it is not possible to bound the number of generators of the sub- 
modules of Sf, that we have to compute with. For instance, for d — 1 and v — 0, choose a positive 

integer k and consider the sub-So-module y^k of Sq generated by the family {~k 3 u 3 )j—Q k- Then 

can not be generated by less than k + 1 elements. Indeed, let (eo, . . . , e„) 6 Sq be a family of 
generators of j$k, and for j > and define a filtration on .Mk by letting F 3 \Jtk = *<$k H u 3 S Q . We 
are going to prove by induction on t 6 {0, . . . , k} that there exists a matrix M t € M nxn (So) such 
that, if we set (e' , . . . , e' n ) = (eo, . . . , e n ) ■ M t then (e' ,...,e' n ) is a family of generators of for 
j < t, e'j = u 3 n h ~ 3 mod F 3+1 Jtk and (e'j)j>t is a family of generators of F l Jik- This is obviously 
true for t = 0. Suppose that it is true for to £ {0, . . . , k}. Let (e , . . . , e' n ) — (e , . . . , e„) ■ M to . As 
the morphism (Y^j—t Sqc^) / F ta+1 ^k — > n k ^ t D\, defined by u to ^ciiU 1 \-t ao is an isomorphism, 
we can suppose if necessary by renumbering the family (eQ that e' to — u to Tr k ~ to mod F ta+1 ^k- 
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Then, by considering linear combinations of the form e'j — Xe' ta+1 for A £ So for j > to, one can 
obtain a matrix M to+ i satisfying the induction hypothesis for to + 1- Finally, we get n > k. 

A second problem comes from the fact that there is no unique way to represent a module by 
a set of generators. For computational purpose, in order to check equality between modules for 
instance, it is important to have a canonical representation, that is a bijective correspondence 
between mathematical objects and data structures. An example of such a canonical representation 
exists for finitely generated modules with coefficients in an Euclidean ring ((Hj): it is the so-called 
Hcrmite Normal Form (HNF). It is given by a triangular matrix (with some extra conditions) that 
can be computed from an initial matrix M by doing operations on column vectors of M . Even 
if Su is not Euclidean, we could have hoped that such representations still exist for free modules. 
Unfortunately, it turns out that it is not the case. Indeed, in general, there does not even exist a 
triangular matrix form for matrices over S v . For instance, for v = 0, take: 

M =( U n ) e A/ 2x2 (5 ) 

\7T UJ 

and assume that M can be written as a product LP where L is lower-triangular and P is invertible. 
Let a and /3 be the diagonal entries of L. Then, a and j3 belong to the maximal ideal of So (since 
the coefficients of M all belong to this ideal) and the product a(3 is equal to a unit times u 2 — it 2 . 
Hence, by multiplying /3 by an invertible element in So if necessary, we can assume that /3 = u ± 7r 
since Sq is a unique factorisation domain. On the other hand, by hypothesis, there exist a, b S Sq 
such that ua + irb — and ira + ub = /?. This equality implies that 7r divides a and therefore that 
/3 = 7ra + ub € uSo + tt 2 So- This is a contradiction. 

In this section, we explain how to get around these problems. First, we recall the notion of 
quasi- isomorphism and study the localisation of the modules with respect to n or u a /tt 13 in order 
to obtain canonical representations well suited for the computation in the category of modules up 
to quasi-isomorphism. Then, we describe a generalisation of an algorithm of Cohen to compute the 
maximal module associated to a given torsion-free S^-module and obtain a bound on the number 
of generators of a maximal 5,,-module. We explain how to combine the different approaches in 
order to obtain a comprehensive algorithmic toolbox for modules over S v . 

3.1 Quasi-isomorphism and maximal modules 

In order to be able to control the number of generators of a S^-module, we are going to compute 
up to finite modules which will be considered as " negligible" . 

Definition 3.1. A finitely generated S v -module is said to be finite if it has finite length. Let ^ 
and be two finitely generated S v -modules, let f : — » be a S v -linear morphism. We say 
that f is a quasi-isomorphism if its kernel and its co-kernel are finite modules. 

Remark 3.2. Since ker/ and coker/ are finitely generated (because S v is a noetherian ring), it 
is easy to check that they have finite length if and only if they are canceled, at the same time, by a 
distinguished element of S v and a power of n. A quasi-isomorphism between torsion-free modules 
is always injective. Indeed, its kernel, being a submodule annihilated by a power of u a /ir^ and tt 
of a torsion free module, is zero. 

Example 3.3. Let .M be the submodule of So generated by (7r 2 ,7ru 3 ). The inclusion C ttSq 
yields an injective morphism whose is annihilated by tt and u 3 . As a consequence ^# is quasi- 
isomorphic to the free module tt.Sq (see figure^). 

We have a canonical representative in a class of quasi-isomorphism which is given by the fol- 
lowing definition. 

Definition 3.4. Let be a torsion-free finitely generated S v -module. We say that together 
with a quasi-isomorphism f : ^ — >• is maximal for ^# if for every JV , torsion-free S v -module, 
and quasi-isomorphism f : „<# — > ,yY , there exists a morphism g : JV — > which makes the 
following diagram commutative: 
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Figure 2: The module ^# is quasi- isomorphic to tt • So- 




(4) 



The morphism g in the definition is unique and is in fact a quasi-isomorphism. Indeed, by the 
commutativity of the diagram, the image of g contains the image of /. Thus, the cokernel of g is 
finite. Moreover, since / is injective, g is injective on Im/', which is cofinite in JV . It follows that 
ker<? is finite and g is a quasi-isomorphism. Moreover, for every x € J\f ', there exists a positive 
integer n such that ir n x is in the image of /'. The image of Tr n x by g is then uniquely defined by 
the commutativity of the diagram (HJ). The uniqueness of g follows. 

A maximal module for if it exists, is unique up to isomorphism. Indeed, if and 
are two maximal modules for then there exist two quasi-isomorphisms g\ : ^# — > and 
(72 : — > and the uniqueness of g in the diagram (jj) implies that g\ o g 2 = Id^» and 
.92 9i = Id^r'. If it exists, we denote the maximal module of ^# by Max(^#). We can rephrase 
the above by saying that if is the maximal module for ^ then there is a quasi-isomorphism 
from into ^f' and any quasi-isomorphism — 1 is an isomorphism. In fact, this condition 
characterises maximal modules: 

Lemma 3.5. Let ^ be a finitely generated torsion free S v -module. Let jH' be a S v -module such 
that there is a quasi-isomorphism f : ,/# — > ./#' . The following assertions are equivalent: 

1. is maximal; 

2. any quasi-isomorphism — > is an isomorphism. 

Proof. We only have to prove that the second property implies that ^#' verifies the universal 
property of maximal modules. For this let JV be a finite type S^-module such that there is a 
quasi-isomorphism /' : ^# — > JV . Let A = / © /' : j$ — > © jV be the diagonal embedding and 
let J(q = -f . It is clear that ^#0 is a finitely generated torsion free ^-module. 



There are canonical injections : — > and i^y : JV — ► ^#o- We claim that 
and % jy are quasi-isomorphisms. To see that, it suffices to show that the induced injection i u 
ifj(' 1 i jv) A : ^# — > ./#o bas a finite cokernel. But 

coker / © coker /' 

coker is/ = 



A(^#) n (coker / © coker/') 



which has finite length being a quotient of coker / © coker/'. 

Next, by hypothesis ijc is in fact an isomorphism so that we have a quasi-isomorphism g = 
i~2, o i y which sits in the following diagram: 



(5) 



It is clear that the lower left triangle of the diagram is commutative and we are done. □ 
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A theorem of Iwasawa [5] asserts that if „4t is a finitely generated module over So, then Max(^) 
exists and is free of finite rank over So- The main object of . 31 is to extend this result to modules 
over S v : we shall provide a constructive proof of the existence of Max(^#) for any finitely generated 
torsion- free module ^# over S v . We will see however that this Max(^#) is not free in general; 
nevertheless we shall provide an upper bound on the number of generators of Max(^#) . 

Lemma 3.6. Let f : j% — > be a quasi- isomorphism between torsion-free finitely generated 
S„ -modules. Suppose that is free then is maximal. 

Proof. We use the criterion of Lemma 13.51 Let JV be a finitely generated S^-module such that 
there is a quasi-isomorphism /' : — > ,yV and we want to show that /' is an isomorphism. As 
is torsion-free, we know that /' is injective. Now, suppose that there exists a non zero element 
in the cokernel of /'. It means that there exists a non zero x € ,JV which is not in the image 
of /'. As /' is a quasi-isomorphism there exists n £ N and A £ S v a distinguished element with 
7r™ ■ x £ Im/' and X-x 6 Im/'. If we set z\ — f 1 (ji n ■ x) and z% = /' 1 (A • x), we have the relation 

\ Zl - tt"z 2 = 0, (6) 

in Let (ej)j 6 / be a basis of ^#' and write Zi — "Ylil l \ e o f° r * = 1, 2. Putting this in we 
obtain that \p,{ — tt 71 ^ and thus ir n \[i{ for j e J since A is a distinguished element of S v . But 
then f'(J2 Mi/ 71 "™^') = ^-/ 7! ' n -f( z i) — x contradicting the fact that x is not in the image of /'. □ 

Remark 3.7. One can rephrase Iwasawa's result in a more abstract way using the category lan- 
guage. Let Mod g be the category of finitely generated S w -modules, that are torsion-free and let 
Mod g (resp. Free ^ ) denote its full subcategory gathering all torsion-free modules (resp. all free 
modules). We also introduce the category Mod g ls , which is by definition the category of finitely 
generated S v -modules up to quasi-isomorphism, i.e. Mod g ls is obtained from Mod g by inverting 
formally quasi-isomorphisms. We have a natural functor Mod g — > Mod g"' , whose restriction to 
Mod g defines a pylonet in the sense of JE/, §L It follows from the results of loc. cit (see Corollary 
1.2.2) that the Max construction is a functor: to a morphism f : .<4t — >• in Mod g , one can 
attach a morphism Max(/) : Max(^) — > Max(^'). We recall briefly the construction o/Max(/). 
Let ./#" be the pushout .Ji' ®j( Max(^), that is the direct sum © Max(^#) divided by ^# 
(embedded diagonally). We have a natural morphism j^C — \ which turns out to be a quasi- 
isomorphism. Hence, there exists a map — > Max(^#') and we finally define Max(^#) to be the 
compositum Max(^#) — > jtft" — > Max(^#') where the first map comes from the natural embedding 
Max(jt) ->• J(' © Max(^). 

If M is a submodule of S% (for some positive integer d) , the following proposition gives a very 
explicit description of Max(^#). 

Proposition 3.8. Write v = f3 /a, with a, (3 relatively prime integers. Let d be a positive integer 
and ^# be a submodule of S„. Then Max(^) exists and 

Max(^) = { x G S[t | 3n e N, n n x e Jt and (tt a //) n • x e Jt }. 

Furthermore the morphism ij( : «/# — > Max(^#) is the natural embedding. 

Proof. Let ./# max be the set of x G S 1 ^ such that there exists some n such that it n x and (u a /tt^) u -x 
belong to j$ . We want to show that Max(^) exists and is equal to ^# max - It is clear that 
j% C Jmax and that the quotient ^# max /^# is canceled by a power of 7r and a power of u a /n 13 
which is a distinguished element. Hence it has finite length, and the inclusion — > ^# max is 
a quasi-isomorphism. Next, suppose that we are given a S^-module together with a quasi- 
isomorphism g : ^# m ax —* ^o- Then there is a quasi-isomorphism ij( : «4( — > that sits in the 
following diagram: 



Jt ► ^max ► fitf 
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Note that g is injective as it is a quasi-isomorphism. Moreover, we know that the cokernel of ij( is 
annihilated by a power of u a /ir^ and a power of ir, which implies that g is surjective. Thus, g is an 
isomorphism and by Lemma l3.5[ Max(^#) exists and Max(^#) = ^# ma x as claimed. The second 
part of the proposition is clear from the above diagram. □ 

It follows directly from Proposition l3.8l that the intersection of two maximal modules is maximal. 
The same is however not true for the sum: in general the S^-module ^# + ^f' is not maximal even 
if and are (take for example M — uSq and Jt' = wSq). This leads us to define the 
new operation + max (which is much more pleasant than the usual sum of modules) on the set of 
maximal submodules of S d as follows: 

J£ + max JC' = Max(^ + J#'). 

We also deduce from Proposition 13.81 that a So-module ^# is free if and only if ^ = ^ max - 
This gives a nice criterion to check if a <So-module is free. It is not true in general for a sub-S 1 ,,- 
module ./# of S d that Max(^#) is free (this will become apparent when we give the general shape 
of a maximal S^-module in ^3.3|) . However, by Lemma 12.51 every S v becomes isomorphic to So 
over a finite extension 9t' = (where w depends on v). Set S' v = S v For all submodule 

j& of S d , we obtain that Max(^ ® Si) is a free submodule of (S' l/ ) d . Denote by Maxg the set of 
maximal sub- S v -modules of S d and by Free^/ the set of free sub-S^-module of (S' 1/ ) d . 

Proposition 3.9. The natural map 

$ : Max^ — > Free^, 

J( i ^ Max(!# ® Sv S' v ) 

is injective. A left inverse of $ is given by jft' H> n S d . Moreover, the image of $ contains 
the subset of Free|f, of free modules which admit a basis (ei)igj where G {Sl) d and = cc ai e^ 
with e- € {S v ) d and a t G N. 

Remark 3.10. Actually, we will prove later (see Lemma \3.18\) that the image of $ is exactly the 
subset o/Free^, verifying the condition of Proposition^^ 

Proof. In order to prove that $ is injective, it is enough to prove that $ has a left inverse. For 
this, let Jt e Maxf^ and let Jt 1 = Max(^# ® Sv S'J) G Free^, . Then it suffices to prove that 
J%2 = H S d is a maximal sub-SV-module of S d . Indeed, as it is clear that ^#2 contains M 
and that the injection ^# — > ^#2 is a quasi-isomorphism since the injection ^# — ¥ is a quasi- 
isomorphism, we remark that by the maximality of it would imply that ^ = ^#2- 

For this let x G S d and suppose that there exists n G N such that 7r" G ^#2 and (u a /ir l3 ) n -x G 
^2- As is maximal and ^#2 C by Proposition 13.81 it means that x G . hence x G -#2- 
Using again Proposition [3TU we deduce that ^2 is maximal. 

Let us now prove the last claim of the proposition. Let G Freeg, which admits a basis 
(ej)j g j where a G (S^) d and e, = nj^ej with e- G and G N. We have to find a sub-S^- 

module of S d such that M®s v S' v is quasi-isomorphic to As M' = ®e.iS' v , it is enough to 
treat the case d = 1. Let < ai be an integer and let be the sub-S^-module of S' v generated by 
w ai . Let A be a positive integer such that ^ + Aj = 7£l Such a A exists because a and ft are 

relatively prime. Let j% be the sub-S^-module of S v generated by 7r and Let = zu~ ai it 
is clear that v„(fj,) = so that is a distinguished element of S 1 ^. Thus, we have ■uj ai ./i G ^ ®s v S' v 
and w" 1 • cc7 Q_Q!1 G ®s u S' v therefore j$ ®s v SI is quasi-isomorphic to . □ 

3.2 An approach based on localisation 

Wc have seen that in a class of quasi-isomorphism of a finite type torsion-free S^-module ^# there 
exists a distinguished element Max(^). In this section, we use this fact in order to represent 
the quasi-isomorphism class of ^ by localizing with respect to u a /tt^ and tt. We thus obtain a 
representation of finite type torsion-free S^-modules amenable to computations. 
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3.2.1 A useful bijection 



We keep our fixed positive integer d. We recall that 

§ = { y diii 1 , at £ K, vk{o-l) + vi bounded below and lim VKidi) + vi = +00) 

is a field containing S v ^ and S u _ u . If J( is a sub-S^-module of <? d , we shall denote by JK^ (resp. 

the sub-S^-module (resp. the sub-S^-module) of S d generated by . For example, if M is 
free over S v with basis (ei, . . . , e/,), then (resp. ^# u ) is also free over S V<E (resp. S v , u ) with the 
same basis. As -M is torsion free, and as S v<w and S v ^ are principal ideal domains, jM v and are 
free. We denote by Maxg^ the set of maximal sub-S^-modules of S d and for A = S V ,S V ^ or S UjU , 
let Free^ denote the set of sub- A- modules of A d , which are free over A. Recall that Max So = Free So . 
Thus, the following lemma provides a useful description of maximal So-modules. 

Lemma 3.11. Let S = So- The natural map 

: Free 5 — > Free^ x Free^ 

is injective. If a pair (A, B) is in the image of ty' , its unique preimage under is given by An B. 

Proof. From the descriptions of elements of 5, S^, S u and S in terms of series, it follows that 
S = SV flS„. If ^ £ Freeg, ^ is isomorphic to S h for h < d and, by applying the preceding remark 
component by component, we get jj( = n^#„. This implies the injectivity of ^' and the given 
formula for its left-inverse. □ 

Using Lemma l3.11[ we can prove: 

Theorem 3.12. The natural map 

* : Max Siy — ► Free^ n x Free^ u 

Jt h-> (^rjJ^u). 

is injective and its image consists of pairs (A, B) such that A and B generate the same $ -vector 
space in S d . If a pair (A, B) satisfies this condition, its unique preimage under ^> is given by Ap\B. 
Furthermore, we have the following equalities: 

^\J(r\J(') = {Jf* n .# u n Jf'^ 
*(^#+ max ^') = (^ + Jf'^, J( u + Jf„) 

for all Jt, Jt' e Max s ^ . 

Proof. Let w in an algebraic closure of K, be such that w a = 7T. Let DV = %K\vj\ and S' u = S v ®fx%K' . 
We know by Lemma f2.5l that SI is isomorphic to IK' [[it]]. Then, the map ^> sits in the following 
commutative diagram: 

* 



Max? > Free? x Free? 



Max(. ® s „ S' v ) 



SL (8) 



Free s / > FrcCg, X FreeJ, 



By Proposition 13.91 the map j$ yisxi^Jt ®s u S' v ) is injective and *]/' is injective by Lemma T3. Ill 
from which we deduce that ^ is injective by the commutativity of ©. 

We want to prove now that if the pair (A, B) belongs to Free s x Free Sy and satisfies the 
condition of the theorem, then = A n B is maximal over S v and 4"(^#) = (A,B). We claim 
that there exists a basis (ei, . . . , eh) of A (over S Vj7T ) such that is included inside the S^-modulc 
generated by the e^'s. Indeed, let us first consider (ei, . . . , eh) a basis of A and denote by the 
5,,-module generated by the e^'s. Now, remark that, by our assumption on the pair (A,B), every 
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element x £ B can be written as a ^-linear combination of the ej's. Taking for n the smallest 
valuation of the coefficients appearing in this writing, we get x £ 7r~"^#^. Moreover, since B is 
finitely generated over S v , u , we can choose a uniform n. Replacing ej by it~ n e' i for all i, we then 
get A = JC!„ and S C JC U . Thus i' = AnBciJniJ| = 

Since S v is a noetherian ring (recall that v is rational) , we find that Jt is finitely generated over 
S v . Furthermore, one can compute Max(^#) using Proposition 13.81 if x is an element of Sf, for 
which there exists n such that ir n x and (u a / n^) n x belong to . // . then x £ A (since it is invertible 
in S v ^) and x £ B (since u a /it? is invertible in S v _ u ). Thus x £ Jt and Max(^#) = Jt, i.e. ^# 
is maximal. 

Let us prove now that ^(Jt) = (A, B). By the same argument as before, we find that there 
exists a positive integer n such that n n ^' C Jt C from what it follows that Jt n = Jt' v = A. 
The method to prove that Jt u = B is analogous: we first show that there exists a basis (ei, . . . , e^) 
of B ov6r Si/ u Emd some elements 5i, . . . , s^, £ Sy such that: 

• all s^s are invertible in S^u, and 

• we have SiCiS v c ^# C X] e i<Sv- 

From these conditions, it follows that Jt u is generated by the ej's over S u and, consequently, that 
Jt u = B. 

It remains to prove the claimed formulas concerning intersections and sums. For the intersection, 
we note that if Jt n Jt' = (Jt v H Jt u ) H n = {Jt* n ..#^) (~1 n ^). Hence, we 
just need to justify that Jt^ D Jt' v and ^„ n Jt' u are free over S v<w and S^m respectively, and 
that they generate the same (^-vector space. The freedom follows from the classification theorem 
of finitely generated modules over principal rings, whereas the second property is a consequence of 
the flatness of $ over S v ^ and S v<w . 

For the sum, we have to justify that (-Jt + m&yi Jt'). n = -^jt+-^ and {Jt + raa ^Jt') u = -^t»+-<^- 
It is clear that (./# + „<#% = + Ji'^ and (•/# + .#')„ = + Jt^. Hence, it is enough to prove 
that, given a finitely generated S^-module N 6 S$, we have Max(iV)^ = and Max(iV) u = iV tt . 
It is obvious by Proposition 13. 81 □ 



Reinterpretation in the language of categories We introduce the "fiber product" category 
Free c^ <8> Frec „ Free g ^ ^ whose objects are triples (A,B,f) where A £ Free s , B £ Free g ^ ^ 
and / : $ ®s v „ A —> S ®s v B is an <?-hnear isomorphism. We have natural functors in both 
directions between Max ^ and Free g ^ <8> Frcc g Free ?i ^ : to an object Jt of Max^ , we associate 
the triple [S v ^ ®s *dt ■> S v ,u ®s Jt , /) where / is the canonical isomorphism, and conversely, to a 
triple {Jt^, ^C, /), we associate the fiber product of the following diagram (which turns out to be 
free of finite rank over S„): 



(9) 



Theorem 13.121 then says that these two functors are equivalences of categories inverse one to 
the other. Actually, this result can be generalized to non-free modules as follows. 

Proposition 3.13. The functor Mod Su — > Mod s ^ g ® Mod g Mod 5 u , Jt i-> {S v ^®sJt 1 S v ^ u ®sJt) 
factors through Mod o' s and the resulting functor 

Mod ^ s -> Mod s, „ ®Mod f Mod o^ u 

is an equivalence of categories. 

Proof. Left to the reader. □ 
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3.2.2 Normal forms for modules over S vw and S v -, 



As S u ^ and S UtU are Euclidean rings there exists a good notion of rank as well as Hermite Normal 
Forms for matrix over these rings. In this section, we state propositions giving the shape of Hermite 
Normal Form together with algorithms with oracles to compute them. We recall that an algorithm 
with oracle is a Turing machine which has access to oracles to store elements of the base ring and 
perform all usual ring operations: test equality, computation of the valuation, addition, opposite, 
multiplication and Euclidean division. We will measure the time complexity of the algorithms by 
counting the number of calls to the oracles. Classically, we then derive some consequences which 
will be used in this paper. For the complexity analysis, we denote by 8 a real number such that 
product of two d x d matrices with coefficient in S v can be done in 0(d e ) ring operations. With 
a naive algorithm, we can take 9 — 3 and with the current best known algorithm of Coppersmith 
and Winograd 0,6* = 2.376. 

Proposition 3.14. Let M = (rriij) £ Mdxd'(S u ,-n), let r be the rank of M . Then, there exists an 
invertible matrix P such that M.P = T with 



T : 







(10) 



y * * J 



where 

• for i = 1 , 



,r, U = u dj + J2iLo 1 °j u ^ w ith v K(bj) + v(j — dj) > 



• for i — 1, . . . , r, Ti^ i = f, and I is a scrictly increasing function from {1, . . . r} to {1, . . . , d} 
such that 1(1) = 1. 

The matrix T is said to be an echelon form of M . Letd-ca^ be the maximal Weierstrass degree of the 
entries of M , an echelon form of M can be computed in 0(d-d' •d max + max(c? e -d', d' e -d) log(2d'/d)) 
ring operations 

If the echelon form moreover satisfies: 

• all entries on the l(i) th -row are elements of K[u] of degree < di. 

then T is unique with these properties and is called the Hermite Normal Form. The Hermite 
Normal form of M can be computed from an echelon form of M at the expense of an additional 
0(r 2 ) ring operations. 

Proposition 3.15. Let M 6 Mdxd'^Sv^u), let r be the rank of M. Then there exists an invertible 
matrix P such that M.P = T and 



T 







\ 







(11) 



where 



• for i — l,...,r, TWjVj = 7r di where I is a strictly increasing function from {l,...r} to 
{1, . . . , d} such that 1(1) = 1. 

The matrix T is said to be an echelon form of M . An echelon form of M can be computed in 
O(d.d') + max(d e • d', d' e ■ d) log(2rf'/d)) ring operations. 
If the echelon form moreover satisfies 



15 



• the entries on the l(i) th -row are representatives modulo ir di . 

then T is unique with these properties and the called the Hermite Normal Form of M . The Hermite 
Normal form of M can be computed at the expense of an additional 0{r 2 ) ring operations. 

Proof. The proof of the previous propositions as well as algorithms to compute the echelon form 
of M with the given complexity is an immediate consequence of Theoreme 3.1] together with 
the fact that S v>w and S^ u are Euclidean rings. Moreover for all x, y £ S v , w one can compute 
the gcd(x, y) in 0(deg w (jj)) ring operations. From its triangle form, one can then compute the 
Hermite Form of M with coefficients in S u _ n at the expense of 0(d ■ r ■ d max ) ring operations. □ 

Remark 3.16. We deduce from this proposition that if M £ Mdxd'(S Ut7r ) is a full rank matrix, 
there exists P such that M ■ P is a matrix of the form \10\) with all coefficients in K[u\. In the 
same way, if M £ Mdxd'(S^ u ) is a full rank matrix then there exists an invertible matrix P such 
that M ■ P has the form f77]) with all entries defined modulo 7 r max { d i>~-> d r} _ 

Let S u ^ioc be S VtU or S u ^. We derive some consequences of the existence of triangle forms and 
Hermite Normal Form for the representation and computation with finitely generated sub-S , t/j ; oc - 
modules of S d loc . We can represent a finitely generated sub-S^oc-module jM of S d loc by a d x d 
matrix M giving d generators of M in the canonical basis of S d loc since every sub- module of S d loc 
has dimension at most d. Keeping the same notations, one can compute the module of syzygies 
of For this it is enough to compute R, a matrix of maximal rank such that M ■ R — which 
can easily be done by computing an echelon form of M. Given a vector ~f £ S~j oc provided by its 
coordinates vector V in the canonical basis, one can check efficiently if f £ jtft by finding a vector 
X such that M.X = V which can also be done with the echelon form of M. 

Let M and M' representing the modules ^# and one can compute a matrix representing 
the module + jM' by computing the echelon form of the matrix (MM 1 ) and taking the d first 
columns. One can compute the intersection of ^# and in the same way by finding R and R' 

such that (MM') ^ J = 0. 

3.2.3 Consequences for algorithmics 

In view of the results of tj3.2.1l and H3.2.21 we shall represent a maximal S^-module ^ living in 
some S d as a pair (A, B) where A (resp. B) is the matrix with coefficients in Sv,-k (resp. in S v ^u) 
in Hermite Normal Form representing S u ^ ®s v (resp. S UtU ®s v 

The second part of Theorem 13.121 tells us that it is very easy to compute intersections and 
"maximal-sums" of S^-modules with this representation. Indeed, we just have to perform the same 
operations on each component, and we have already explained in H3.2.2I how to do it efficiently. 
As the Hermite Normal Form is unique, it is also very easy to check the equality of two maximal 
sub-S^-modules of S d . Using only the echelon form of the matrices A and B it is also possible to 
test membership. 

Even better, this representation is also very convenient for many other operations we would 
like to perform on 5,,-modules. Below we detail three of them. First, let ^# C S d be a maximal 
S^-module. By definition, the saturation of ^# in S d is the module 

^ sa t = { X £ S d | 3n £ N, n n x £ Ji). 

It follows from Proposition ^. 81 that ^# S at is maximal over S v , and we would like to compute it. For 
that, working with our representation, we need to compute (^sat)n- and (./# sat ) u . But, we have 
(-^sat)*- = J(-k and 

(^sat)« ={16^, I 3n £ N, 7r n x £ JZ U ). 

The computation of (^sat^ is then for free, whereas the computation of (^# S at)« can be achieved 
using Smith forms, which is here quite efficient due to the fact that S V}U is a discrete valuation 
ring. An important special case is when ^# has rank d over S v . Then (^# sa t)« is always equal 
to S d u . Thus, in this case, if ^# is represented by the pair of matrices (A, B), then ^# sa t is just 
represented by the pair (A, I) where I is the identity matrix. 
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More generally, one can consider the following situation. Let j$ g Maxg and jtf,' G Maxg 
We want to compute ^#n^#', which is a maximal module over S v . As before, we need to determine 
n ^C')^ and n ^#') u and one can check that: 

(M n M\ = Jin n JC V 

(JtC\Jt') u = Jt u C\Jt' u . 

Note that, here, ^i' u is vector space over S . As before, the intersection ^#„n^#^ can be computed 
using Smith forms and, if has rank d over S v>n , we just have j$' n = <? d and so (^V\^') u = 

The third example we would like to present is obtained from the previous one by inverting the 
roles of S^tt and S ViU : we take € Free^ and ./#' € Free^ and we want to compute ^# n 
We then have (j$ n = n and n ./#')„ = ^ u fl Here a new difficulty occurs: 

is a (^-vector space and so, in previous formulas, it appears an intersection between a free 
module over S v ^ and a <?-vector space. Again, one can compute this Smith form. However, it 
is not so efficient as before since S v ^ is just an Euclidean ring, and not a discrete valuation ring. 
Anyway, it remains true that, in the case where has full rank, then = S d . So, in this case, 
n is just equal to and the computation of (Ji( n becomes very easy. 

3.2.4 Further localisations 

We remark that the matrix appearing in Proposition [3715] has coefficients in S VM which is a discrete 
valuation ring while the matrix of Proposition 13 . 141 has coefficients in S v ^ which is only Euclidean. 
For certain applications, it can be more convenient to compute with elements in a discrete valuation 
ring; for instance, the computation of the Smith Normal Form can be made faster in a discrete 
valuation ring. 

It is actually possible to work only over discrete valuation rings by localising further. More 
precisely, for any element a € K (where K is an algebraic closure K of K) with valuation > v, we 
have a canonical injective morphism S VfW — > K[[u — a\] which maps a series to its Taylor expansion 
at a. Hence, if is a sub-S^-module of S d n , one can consider K[[u - a}} C 

K[[u — a]] d for all element a as before. Moreover, if ^# p has maximal rank, all ^# PjQ 's are trivial 
(i.e. equal to K[[u — a]]) expect a finite number of them (which are those for which a is a root of 
one of the V s of Proposition 13 . 14|) . In addition, the map: 

3 :Mod^ — ► n^ M "4[[ u -«]] 

Jltp l-> (J( v ,a)a 

is injective and commutes with sums and intersections. Hence, one can substitute to ^# p , the 
(finite) family consisting of all non trivial ^# p>a 's. This way, we just have to work with modules 
defined over discrete valuation rings. 

Note finally that there exist algorithms to compute one representation from the other. Indeed, 
remark first that computing the image of ^# p by 5 is trivial if is represented by a matrix of 
generators: it is enough to map all coefficients of this matrix to all K[[u — a]]'s. Going in the other 
direction is more subtle but is explained In [3J, §2.3. 

3.3 A generalisation of Iwasawa's theorem and applications 

The aim of this subsection is to present an algorithm with oracle to compute the maximal module 
associated to a S^-module. Moreover, as a byproduct of our study, we will derive an upper bound 
on the number of generators of a maximal sub-S^-module of S™ ■ 

The idea of our construction (inspired by an algorithm of Cohen) is to consider the matrix of 
relations of a module and to perform elementary operations preserving quasi-isomorphisms to put 
this matrix in a certain form. In order to do so, we first need a way to compute the matrix of 
relations of a module or at least a certain approximation of it. Let ^ be a torsion-free finitely 
generated S^-module and let (ei, . . . , e&) G ^ k be a family of generators of We denote by 
the module of relations of (ei, . . . , e^) that is the set of (Ai, . . . , Afc) £ S k such that J2i=i ^i e i = 0- 
Let r be the rank of ^# ®s„ S Ut „. From the exact sequence 

-> M ® Sv S„,* -»• S*,, -> Jt ®s v S V<K -> 0, (12) 
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we deduce that Si ®s„ S v ^ is a free module over S V<7V of rank I = k — r. Let . . . , /e) be a 
basis of ^ ® s „ S^tt and set ^" = ©1^(5^ ■ /jflSj). Apparently, ^' is a sub-S^-module of ^ 
which is free of rank £. Indeed, if denotes the smallest integer such that 7r ni ■ fa G S*, then 
the family (7r ni • /j) is a basis of S?'. Moreover, we have the inclusion Si' D ir N S? for a certain iV 
since ^" ®g u = ^ ®s„ S Vt7r . Now, from the knowledge of the matrix M € Mdxk(Su) whose 
column vectors are the coordinates of e, in the canonical basis of S„, we can compute a matrix 
R' £ Mkxi{S u ) of generators of Si' using the algorithms of M3.2.2I We have by definition M.R' = 0. 
Of course in the above construction, we can replace, mutatis mutandis the localisation with respect 
to 7r by the localisation with respect to u a /n 13 . 



3.3.1 An algorithm to compute the maximal module 

We start with a couple of matrices M = (rriij) € Mdxk(S v ) and R = (r^j) 6 Mf. x i(S v ) representing 
the generators of ^# embedded in Sf, and a sub-module of Si containing -k Si for a certain N. 
We are going to prove by induction that we can put R in triangular form by using elementary 
operations on the rows of R and the columns of M which preserve ^# up to quasi-isomorphism. 
We suppose that for a positive integer io there is a strictly increasing function t : [l,io] N* such 
that 



• for all i = 1, . . . , io 

• for all i = 1, . . . , io, for all j > t(i), rij 
The matrix R has the following shape: 

/ 



1, for j > i, and t(i) < m < t(i + 1), r JiT 
0. 



= ; 



R 



ri,tm 



r i ,t(i ) 



(13) 



v *-*; 

where the blanks represent entries. 

We set t(io + 1) to be the first integer t such that t(io) < t < I and there exists a j > io + 1 
with Tj t t 7^ 0. If no such integer exists then we have finished. In order to describe operations on 
rows (resp. columns) of a matrix T of dimension k x I it is convenient to denote the row vectors of 
T (resp. the column vectors of T) by Li{T) for i = 1, . . . ,k (resp. Ci(T) for i = l,. ..,£). We say 
that the condition Cond(i) on R is satisfied if there exist two different indices jo,ji G {1, . . . , k} 
such that r JoAi) -r jlAi) ^ 0, v v (r J(ut{i) ) < v v (rj ut (i)) and deg w (r joit ^) < Aeg w {r Jut{€} ). We apply 
the algorithm ColumnReduction (see Algorithm |3]) on R, M, io + 1, t(io + 1). 

It is clear that the matrix M returned by Algorithm [3] represents the same module since 
it modifies M by performing elementary operations on the columns. Moreover, the algorithm 
preserves the relation M ■ R — 0. The effect of the operation of Step [5] of Algorithm [3] on the entry 
rji tu) °f R is either 

• replace it by 0, 

• or it decreases strictly its Weierstrass degree and it increases its Gauss valuation. 

Hence, it is easily seen that after a finite number of loops the conditions Cond(£(io + 1)) will 
no longer be satisfied on R. It may happen that there is only one nonzero entry on the t(io + l)*' 1 
column of R and in this case, we are basically done: by permuting the rows of R we can suppose 
that the non zero entry is ^i +i.t(i +i)- Next, we remark that the vector v of ^ whose coordinates 
in the canonical basis of S„ is given by the (io + l) th column of M verifies r io+1 ,t(j +i) -v = which 
means that v = and we can set r^+ij = for j > t(io + 1). 

If there are several nonzero entries on the t(io + l) th column of R and the condition Cond(t(io + 
1)) is not satisfied on R, we let jo be such that Vv(rj ,t(i +i)) = m ^ n -i<j<k{Vi/(fj v t(i +V))}- Note 
that we have ^v(fj ,t(i +i)) < v v( r j,t(i +i)) f° r 3 ^ io because on the contrary, the condition 
Cond(t(io + 1)) would be satisfied on R. By multiplying the t(io + l) th column of R by an element 
of S V} n with valuation —Vu(fj 0! t(i +V))i we can moreover suppose that v„(rj g t u 0+ -n) = 0. Let 
S = rmn j¥:jo (v u (r jAio+1) )). 
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Algorithm 3: ColumnReduction (preliminary version) 



input 

• Me M dy . k {S„) 

• Re M k xi{S v ) in the form ((T5J}, 

• i,t(i) e N 

output : R,M such that M ■ R = and i? does not satisfy condition Cond(t(i)) 

1 while Cond(i(i)) is satisfied do 

Pick up j ,j'i € {1, . . . ,fc} such that r\, M(l) • r jl)t(i ) ^ 0, ^(r,^)) < w„(r ilit(l0+1) ) and 
deg w (r joMt) ) < deg w (r Jut{l) ); 
(q,r) <- EuclideanDivision^^;),^^)); 
C ]0 (M)^C JO (M)+qC n (M); ' 
L h (R) <- L h {R) - qL jo (R); 

6 return M, R: 



The case v = First, we suppose that v = from which we deduce that 5 is a positive integer. 
Denote by ei,...,ek the generators of ^ represented by the column vectors of the matrix M. 
Denote by the module generated by (e'j)j—x...k with e'j = ej for j ^ j and e^ = \e.j a - The 
identity of S 1 ^ induces an inclusion / : ^# — \ It is clear that the cokernel of / is annihilated 
by 7r. Moreover, we have 

R r ( 14 ) 

As the right hand side of (fPfj) is in ^ since rj - t( ^° +1) g S v , the cokernel of / is also annihilated by 
Tj 0t tu +x) which is a distinguished element of S v . We conclude that / is a quasi-isomorphism. 

We denote by 0\ (j) the operation on the couple of matrices (M, R) which consists in multiplying 
by i the (j) th column of M and multiplying by ir the (j) th row of R. Keeping the hypothesis and 
notations of the preceding paragraph, it is clear that if (M, R) represents the module ^ and its 
relations, then the matrices resulting of the operation of Oi(jo) represents the module which 
is quasi-isomorphic to ./#. By repeating operations of the form 0%(j) a finite number of time, we 
can suppose that S = 0. But it means that the condition Cond(i(io + 1)) is not satisfied on R and 
we can call again Algorithm [3] 

We thus obtain the algorithm ColumnReduction (final version) which takes a relation matrix 
of the form (|13[) for io and returns a relation matrix of the same form for io + 1. The algorithm 
MatrixReduction (final version), Algorithm El uses ColumnReduction in order to compute a new 
set of generators of a module quasi-isomorphic to ^# the relation matrix of which has a triangular 
form. 

The general case We reduce the general case to the case v = 0, by using Lemma 12.51 Let 
w in an algebraic closure of K be such that w a — ir. Let 9t' = Stfnj], S' v — S v <E)fn 9V and 
= ^ ®s v S' v . The valuation on (resp. the Gauss valuation on S^) extends uniquely to 
9t' (resp. to S' v ). We have v v (m) = 1/a. The algorithm for the general case is exactly the same 
as for the case v = up to the point when Cond(i(io + 1)) is not satisfied. By multiplying the 
t(i + l) th column of R by w~ Vl, ^ r ' : >o. t ( i a+^'>' a ; we can moreover suppose that Vv(rj Mi +i)) = 0- Let 

With this setting, we can define a quasi-isomorphism in the same manner as before. Namely, 
let ei, . . . , efc be the generators of ^M' as a sub-module of S' v represented by the column vectors 
of the matrix M. Denote by ^[ the module generated by (e'j)j=i...k where e'j — ej for j ^ j and 
e' o = ^ej . Then the natural injection Jl' — > JK[ is a quasi-isomorphism. We denote by O2O, S) 
the operation on the couple of matrices (M, R) with coefficients in S' v which consists in multiplying 
by ^7 the (j) th column of M and multiplying by w s the (j) th row of R. With the hypothesis and 
notations of this paragraph (i.e. M has the form (fT3"|) ). if (M, R) represents the module jtf' and its 



19 



Algorithm 4: MatrixReduction for the case v = 



input 

• ReM kxe (S u ), 

• Me M dxk (S u ) such that M ■ R = 0. 

output : R e Mkxi(Sl), M G Mdxk(S' v ) such that M • i? = and R is a triangular matrix, 
l i <- 0; 



2 t(i ) <- 1; 

3 while i < k do 



i(*o) ^- min{i|t > t(io) andBj > 0, withr^t ^ }; 

«o <- «o + l; 

M, i? <— ColumnReduction(M, _R, iq, f («o))j 
for j <— t(i ) + 1 to £ do 

L r 'o,j <- 



Algorithm 5: ColumnReduction (final version) for v = 



input 

• M G M dxk (S v ), 

• Re M k xt{S v ) in the form ([15]), 

• i, G N the position of the last non zero "diagonal" entry of R. 

output : R,M such that M.R — and R is triangular up to the i + 1 row. 

1 while 3j ,ji such that j ± ji and r joMi) ■ r jut{i) ^ do 

2 while Cond(t(i)) is satisfied do 

3 Pick up j , ji G {1, . . . , k} such that r joit (i) ■ r jut {i) ^ 0, v u (r j(ut{i) ) < v u {r n , t{i) ) and 
deg w (r JOjt(i) ) < deg w (r Jum ); 
(q,r) <- EuclideanDivision(r i0it(i ),r iltW ); 
C j0 {M)^C jo {M)+qC h {M); 
L h (R) <- L h {R) - qL jo (R); 

Let jo be such that deg w (r j0}t{i )) = maxi< J </ s {deg w ,(r i)t ( i ))}; 

C j0 (M) <- £C j0 (M); 
L jo (R)^7r s L jo (R); 

n return M, i?; 



4 
5 
6 

7 
8 
9 
10 
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relations, then the matrices (M',R') resulting of the operation of 02(30, S) represents the module 
^[ which have been shown to be quasi- isomorphic to (as a ^-module). Moreover, R! verifies 
the condition Cond(t(io + 1)). 

The matrix M' (resp. i?'), resulting of the operation O2G, S) is made of column (resp. row) 
vectors with coefficients in S v multiplied by w s for a certain S £ —Z. An important claim is that 
this structure will be kept intact in the course of the computations involving all the elementary 
operations introduced up to now. In fact, these operations on the rows of R are: 

• multiplication of a row by a w a , for a an integer ; 

• permutation of the rows ; 

• for jo, ji £ {1, . . . , k}, replacing Lj 1 (R) by Lj 1 (R) — q'Lj (R) where q 1 is the quotient of w" -y 
by w a ° ■ x for x, y £ S v and ao, a\ £ N. 

It is clear that the two first operations does not change the structure of R and the same thing is 
true for the last operation. Indeed, let q £ S v ^ and r £ S v ^ n K[u] with deg(r) < deg w (x), be 
such that y — q ■ x + r, then for ao, ol\ £ N, we have w a ° ■ y — zu a °~ ai q ■ m ai x + zu a °r so that we 
have q' = w a °~ ai q with q £ S v . 

In order to prove formality this claim and take advantage of it to carry out all the computations 
in the smaller S v coefficient ring, we represent the couple of matrices (M',R') with coefficients in 
SI by a triple (M, R, L) where M, R are matrices with coefficients in S v and L = \a\, . . . , a*] is a 
list of integers such that for i = 1, . . . , k, Ci{M') = rof C,-(M) and L^R 1 ) = w^L^R). We say 
that the condition Cond'(z) on R is satisfied if there exists two different jo,ji £ {1> • • • , k} such that 

r j ,t(i) • r n,m ^ °> M r h,m) + - v »( r h,m) + and de sw( r j«A^ - de &w( r h,m^ With 

these notations, we can write the final version of the MatrixReduction algorithm (see Algorithm 
E|) which encode the matrices M',R' with coefficients in SI with a couple M,R of matrices with 
coefficients in S v and a list of integers. 

Example 3.17. We illustrate the operation of the algorithm on the module of example \3.3{ Recall 
that ^ is the submodule of Sq generated by (7r 2 ,7ru 3 ). It is represented in the canonical basis of 
So by the matrices M of generators and R of relation : 

M = (vr 2 7tu 3 ) , R = 

It is clear that Cond(l) is not verified on R since there is no division possible between its entries. 
As a consequence, we apply operation Oi(l) on the couple (M,R) to obtain: 

Now, we have iru 3 = —u 3 ■ ir and by applying on M (resp. R) an elementary operation on the 
columns ( resp. rows ), we get finally : 



M = (tt Q),R = 




An we deduce that the maximal module associate to ^ is tt.Sq- 
3.3.2 Computation of Max(^f) 

Let Mi,Ri,Li = MatrixReduction(Af, R, L = [0, . . . , 0]). Let L x = [ft, . . . , p k ]. We denote by J([ 
the sub-S^-module of (S' v ) d generated by the vectors given in the canonical basis of {S' v ) d by the 
column vectors ■ Ci{M\) for i £ {1, . . . , k} such that Li(Ri) is the zero vector. 

Lemma 3.18. We have Jt[ = Max(^# ® Sv S'„). 
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Algorithm 6: MatrixReduction 



input 

• ReM kxe (S v ), 

• Me M dxk (S u ) such that M ■ R = 0. 

output : R e M kxi (Sl), M G M dxk {S' v ), L such that M ■ R = and i? is a triangular 
matrix. 

1 i <- 0; 

2 t(i ) <- 1; 

3 L [0, . . . , 0]; 

4 while i < k do 

«o «- «o + 1; 

t(io) <r- min{t|t > t(io)&nd3j > 0, withr^ ^ }; 
while 3.7o,. ?i sucft that j ^ ji and r joit ( io ) ■ r jut(io) ^ do 
while Cond'(t(i )) is satisfied do 

Pick up jo, ji G {1, • • • , k} such that r jojt(io) • r jl>t(io) ^ 0, 

««'( r jo,t(<o)) + ^ < ^(^i,t(io)) + ^ and de gw( r jo,t(io)) < de Ew( r nA*o)Y> 
if ^(^ ,t(io)) > u "( r ji,t(»o)) then 

So <- KOW(io)) -v v {r jlAio) )\, 
L n {R)^^L n (R)- 

L[ji] <r- L{j{\ +a ■ S ; 

(q,r) <- EuclideanDivision(r joi i(j ),r Jlit(io )); 
C jo {M)^C jo {M) + qC h {M)- 
L h {R) «- L h {R) - qL jo (R); 

Let j be such that deg w (r loMlo) ) = maxi<j< fc {deg w ,(r J - )t ( io ))}; 
S <- rmajjij {v v {r jit{ i o) )) - v v {r jott{io) ); 
C j0 (M) <- ^rC j0 (M); 
L j0 (R)^nl s ±L j0 (R); 
L\jo]<-L\j ]+5-[5\; 
for j «- t(i ) + 1 to £ do 
L r »»J <~ 
Let jo G {1, • • • , k} be such that r j0yt(ia) ^ 0; 
(C j0 (M),C i0 (M)) <- (C io (M),C jo (M)); 
(L jo (i?),L io (i?))^(L io (i?),L jo (i?)); 
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Proof. Let = ^ ®s u S' v and let ^#1 be the sub-S^-module of {S' v ) d generated by the column 
vectors of M\. It is clear that ,M\ = since for i G {l,...,fc} such that Li(Ri) is not the 
zero vector, we have Cj(Mi) = (because ^\ is torsion free). As is obtained from jM' by a 
sequence of quasi-isomorphisms, it means that there exists a quasi-isomorphism q' : jtfC — > If 
we prove that is a free ^-module, we are done by Lemma |3~61 

Consider the exact sequence — »■ 2% — »■ — >■ — »■ associated to the family (ei, . . . , e&) of 
generators of As S' v is flat over £„, and as ffl' ®s u S' v [l/w] =3% ®s v S' v \\/vj\ by definition of 
ffl' , we have an exact sequence 

-> # <g)s„ S^l/ro] -> (5; fe )[l/n7] -»• j£'[l/w] -¥ (15) 

defined by the generators (ei,...,ejt) of It is clear that at each step, the algorithm 

ReduceMatrix describes an exact sequence of the form (|15p for a different map (S' v k )[l / 'vo\ — > 
^'[1/ro] since it preserves the relation MR = 0. From this and the definition of M[, we deduce 
that if &\ is the module of relations of J([ then 9i\ [1 /w] = from which we deduce that 3&\ = 
and we are done. □ 

Remark 3.19. As a byproduct of the preceding proof we see that the vectors given in the canonical 
basis of (S' u ) by the column vectors vj^ i ■ Ci{M\) for i G {1, . . . , k} such that Li{R\) is the zero 
vector form a basis of . 

Corollary 3.20. Let Jl 2 = Jt{ n S*. Then, Jt 2 = Max(^). 

Proof. The corollary is an immediate consequence of Proposition 13.91 and Lemma 13.181 □ 



3.3.3 Computation with 5„-modules 

Proposition 13.91 and Lemma 13.181 establish a one-to-one correspondence $ : Max^ — > Freeg/ , 
defined by >— > Max(^# <E>s v S' v ). Moreover, the image of $ is exactly the set of free sub- 5^,- 
modules of S' v d which admit a basis (e;)i e / where e, G {S' v ) d and e, = , cu oli e' i with e' i G {S v ) d and 
< on < a. We have seen that a j$ G <I>(Maxg ) can be represented by a couple (M,L) where 
M G Mdxk(Su) and L is a list of positive integers < a. 

From the data of a matrix representing an element of ^ G Max^ the algorithm MatrixReduc- 
tion computes the couple (M, L) representing $(^). Moreover, if G $(Max^ ), the Algorithm 
[7] allows to recover $ _1 (^#'). We see that we can easily go back and forth between the different 
representations. For most of the applications however, it is convenient to represent an element of 

G Maxg by a couple (M, L). Indeed, we have the lemma: 

Lemma 3.21. Let Jt\,Jt 2 G Maxf^, then 

$(4nJ 2 ) = *(-*r x ) n $(-#2), 

+ max J^2) = $pfl) +max *(-^2)- 

Proof. For the first claim, we have $- 1 ($(^i) n $(^# 2 )) = n $(^ 2 ) n5^ = ($(^#l) n 

Si) n ($(^ 2 ) n s£) = J( x n -# 2 . 

Next, we prove the second claim. We have the following diagram of quasi-isomorphisms: 




(16) 



Max(^#i ® s „ S' v ) + Max(^# 2 <®s„ S^,) 



Thus, we have Max(Max(^i + Jt 2 ) ®s„ S' v ) = Max((^#i + ^# 2 ) ® Sv S' v ) = Max(Max(^i ® Sv 
S' v ) + Max(^ 2 (g)s v S' u )) which is exactly the desired result. □ 

Let j#\,^2 G ^(Maxg ) be represented respectively by the couples (M\,L\) and (M 2 ,L 2 ). 
Then, by Lemma 13.211 one can represent the sum + max -^i by applying the algorithm Ma- 
trixReduction on the couple ((MiM 2 ),ii + L 2 ) (where L\ + L 2 is the concatenation of the lists 
L\ and L 2 ). The representation as a couple (M,L) is however not well suited to the computation 
of the intersection of modules, since it implies the computation of the kernel of a matrix with 
coefficient in S v which is not Euclidean. 
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3.3.4 The generators of a maximal module 

In order to have a complete algorithm (with oracles) to compute Max(.,#), it remains to explain how 
to recover Mi = ^[C\S^ from the knowledge of Ji[ (see §3.3.2l for the definition of We would 
like also to obtain a bound on the number of generators of ^#2- By the construction of .4ff{, there 
exists a basis (ei, ...,et)eS^ and Si G N for i = 1, . . . , k, such that M[ = 0,_ 1 S' v .w &i ei. Then, 
we have .-#2 — ©^(^.tu^ nft,)-^. Hence, it is enough to explain how to compute ^{nS^ when 
has dimension 1. In this case, is generated by an element of the form -\ ■ y where y G S v 
and by definition, we want to find generators for the ft,-module {x G S v \v v (x) > v u (-\ ■ y)}. We 
are reduced to the problem of finding generators of the S^-module ,jV = {x G S v \v v (x) > — 5/a}. 

Lemma 3.22. Let 5 G {0, . . . , a — 1}. We define inductively a sequence of couple of integers 
(ai,/3i) by setting ao = 0, ft) = 0. Then for i > 0, while ft_i + a^\v > — — , we let (cti,ft) be the 
unique couple of integers such that 

• ft + OLiV > 

• for all (x, y) 7^ (ai, ft) G Z 2 such that < x < ai and y+xv > ——, we have Pi + aiV < y+xv, 

• ai is the smallest integer strictly greater than a^_i such that there exists an integer ft with 
(«i,ft) satisfying the two conditions above. 

The family (7r ,3i • u ai ) has cardinality bounded by a and is a system of generators of the S u - 
module jV — {x G S y \v y (x) > — 6/ a}. 

Proof. First, it is clear by definition that all the ■ u ai are elements of ,JV . Moreover, it is clear 
that ai is bounded by —S//3 mod a. 

Denote by jVq the sub-ft.-module of generated by the family (71-^ ■ u ai ). Let x G jV , we 
prove inductively on deg w {x) that x is in JV§. If deg w (x) = then v v {x) > so that x — x T with 
x G S v . Suppose that d = deg w {x) > 0. As v„(x) > —5/a, by applying Corollary |2.9l we can write 
x = q ■ h, with q G S v invertible and h G K[u] is a degree d polynomial such that v v (h) > —5/a 
and deg w (h) — d. We have to show that h is in jVq. Let io be the greatest index such that ai < d. 
Then by construction of the family (o^, ft), we have w 1/ (7r /3i o • u ai o ) < v u (h). Indeed, if t is the term 
of h of degree d then t G jV and if we write t = 7r M • u x , we have by construction ft + ct; v < fi + xv. 
Thus we can write h = q\ ■ 7r^o • u ai o + r where q\ G S v , deg w {r) < ai and v v {r) > —5/a. We 
can then apply the induction hypothesis on r to conclude. □ 

From the above lemma, one can easily deduce an algorithm to compute the generators of 
JV — {x G S u \v v {x) > — 5/a} as well as an upper bound on the number of generators. In order 
to find the ai we just run over all the values between 1 and —5/(3 mod a and check for each of 
them if it satisfies the conditions of Lemma 13.221 Nevertheless this algorithm is inefficient and 
the obtained bound is far from tight. In the following, we explain how to obtain a tight bound 
as well as an efficient algorithm to compute a family of generators of JV by using the theory of 
continued fractions. In order to set up the notations, we briefly recall the results from this theory 
that we need (see [9]). For ao, . . . ,a n integers, the notation [ao;a±,..., a n ) refers to the value of 
the continued fraction 

1 

a H . 




We take the convention that a n 7^ 1 in [ao; a±, . . . , a n ] so that every rational number can be 
written uniquely as a finite continued fraction. Let r = [ao; 01, . . . , a n \. We let pq = ao, 170 = L 
pi = a ai + 1, qi = a\ and define inductively p k = a fe p fc _i + Pk-2, qk = a>kqk-i + qk-2- The 
fractions p k /qk are called the k th convergent of the continued fraction [ao; a\, . . . , a n \. We have the 
properties: 

• the integers pk and qk are relatively prime (see Th. 2]); 

• Pk/qk = [ao; ai, . . . ,Ofc]. 
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Definition 3.23. Let r be a real number, and let 7 be a positive integer. We say that a fraction 
§ (b>~i) is a best approximation (resp. a positive best approximation ) of r relatively to 7 if for 
all integers c, d such that 7 < d < b and c/d ^ a/b (resp. such that < d < b, dr — c > and 
c/d / a/b), we have \dr — c\ > \br — a\ (resp. dr — c > br — a > 0). We say simply that | is a 
best approximation (resp. a positive best approximation^ of r if £ is a best approximation (resp. 
a positive best approximation) relatively to 1. 

Remark 3.24. Our definition of best approximation corresponds to what is often called in the 
literature best approximation of second kind (see \9j). 

Everything we need about continued fractions is contained in the following theorem (see [9j Th. 
15 and Th. 16]). 

Theorem 3.25. Let x = [do; 01, . . . , a n ]. 

1. Every convergent Pk/qk is a best approximation of x. 

2. Reciprocally, every best approximation of x is a convergent, the only exceptions being the 
cases x = a Q + k, with k G [1/2, 1[, |a = Si. 

Moreover, for i = 0, . . . , n — 1, x — 2i > for i even and x — ^ < for i odd. 

Let r be a real number and b an integer. In the following, it is convenient to denote by min(r, b) 
(resp. min + (r, b)) the integer a such that \b ■ r — a\ — min{|6 -r — k\,k G Z} (resp. such that 
b ■ r — a = min{fe • r — k,k G Z with b ■ r — k > 0}). Then, for r a real number and b a positive 
integer, we let {b} r = b ■ r — min(r, b) and {b}£ = b • r — min + (r, b). 

Example 3.26. Let r — 0.9 and 6 = 2. Then we have min(r, b) — 2, min + (r, b) — 1, {&},. = -0.2 
and {b}+ = 0.8. 

We need the following lemma: 

Lemma 3.27. We have: 

• for all j G {0, . . . ,n}, {q 3 } x > if j is even, {q 3 } x < if j is odd; 

• for j G {1, ... ,n — 2} for all £ integer such that < C < a j+2, C ' {qj+i}x + {qj}x has the 
same sign has {qj} x - 

Moreover for all j G {1, . . . , n — 2} and all £ integer such that < C < a-j+2, 

{C • Qj+i + qj}x = C ■ {q 3 +i}x + {q 3 }x- 

Proof. The fact that {qj} x > if j is even, {qj} x < if j is odd is an immediate consequence of 
Theorem E25l 

If C = 0, there is nothing to prove. We suppose for instance that {qj} x > and {qj+i} x < 
(the other case can be treated in a similar manner). Suppose that for < C < we have 

!'/..! r • . •!'/.,• : K, • 0. (17) 

Let £ be the smallest verifying (|17|) . then £ > 2 since we have by definition of a best approximation 

\{qj}x\ > \{q ]+ i} x \. Then,aa{qj} !C + (C-l)-{qj + i} x >0,weha.ve\{qj} x + C-{qj+i}x\ < \{qj+i}x\ 
which is a contradiction with the fact that there is no best approximation of x the denominator of 
which is between qj+\ and qj+2 = o-n+2qj+i + qj > C ' Qj+l + Qj- 

With our hypothesis, for all integer £ such that < ( < aj + 2, we have {qj} x > {qj} x +(-{qj+i} x . 
Thus we have we have {qj} x > C(qj+i-x— min(a;, qj+i))+qj-x— min(x, qj) > 0, so that 1/2 > (C<Zj+i + 
qj) ■ x — Qm\rt{x, qj+i) — min(a;, qj) > (remember that as j > 1, {qj} x < 1/2). As a consequence, 
( min(x, q ]+1 ) +mm(x, qj) = min(x, (q j+ i + qj) thus {C • q j+1 + qj} x = £ ■ {qj+i} x + {qj}x- □ 

For x = [ao; a±, . . . , a n ] G Q and 7 a positive integer, we would like to be able to obtain 
the list of positive best approximations of x relatively to 7. The lemma tells us that not only 
the convergents P2i/q2i for i G {0, . . . , |/i/2j} are positive best approximations of x but also the 
min + (x, q 2 i + Mfe+iVfe + W2i+i) for i G {0, . . . , [(n — 2)/2j } and \i integer such that 1 < fi < 
&2i+2- The following proposition states that these are all the positive best approximations of x and 
gives a generalisation for the case of a positive 7. 
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Proposition 3.28. Let x = a/b where a, b are relatively prime integers. Write x = [do; ax, ... , a n ] 
and denote by Pk/qk the sequence of convergents associated to the continued fraction [ao; ai, ■ • ■ , a n ]. 
Let 7 < b be a positive integer. Let 7 < d < b be an integer such that min ( x > d ) j s a positive best 
approximation of x relatively to 7. Let i be the biggest index such that d — q 2 i+i > 7 ond let A be 
the biggest integer such that d — <?2i+i — A • <?2i+2 > 7- Tften 

J) mln ( x -. d -<i2z+i->'-q2z+2) j s positive best approximation of x relatively to 7. 

2) If e is such that d— q 2 i+\ — \-q 2 i+ 2 < e < d then min + (x, e)/e is no£ a positive best approximation 
of x relatively to 7. 

Moreover, we have 

{d - q 2i+ i - A • q 2 i+ 2 }t - { d )t = x ' {?2i+2}x - tei+i}* > 0. (18) 

Proof. Let i and A be denned as in the statement. We remark that we have A < 02^+3. Indeed, by 
hypothesis d - q 2i+ i - A • q 2l+2 > 7, but we have q 2l+ 3 = 021+3 ■ ?2i+2 + ?2i+i and we know that 
d - <72i+3 < 7- For < ( < a 2 i+3 an integer, let /Lt(£) = q 2i+ i + C,.q 2 i +2 , h = d- p{\). 
First, we prove that 

{d}i-MOh = {d-»(0}t, (19) 
if < C < a2i+3. Using Lemma 13.271 we obtain 

< mm(x,n(()) -fi(() -x < 1. (20) 

As < d ■ x — min + (a:,c?) < 1, we have < (d — p{()) • x — min + (a:,c?) + min(x, p(()) < 2. We 
have to prove that (d — /i(£)) • x — min + (a;, d) + min(a;, p(()) < 1. Suppose, on the contrary, that 
(d — /i(C)) • x — min + (:r, d) + min(;r, p(()) > 1, then because of (|20|). we have: 

< (d - fi(C)) ■ x - min + (a;, d) + min(a;, /«(C)) - 1 < d ■ X - min + (x, d). (21) 

If ( < A this is a contradiction with the hypothesis that min ( x ' d *> [ s a positive best approximation 
of x relatively to 7. If ( > A then (d — ^i(C)) • x — min + (a;, d) + min(a;, fJ,(Q) < (d — /u(A)) • x — 
mm + (x,d) + min(x, ^i(A)) because {p(C)}t > {m(^)}^ by Lemma 13.271 Next, we remark that 
(d — p(X)) • x — min + (x,d) + min(x,p(X)) < 1 by what we have just proved, so that we have 
(d — p(()) • x — min + (x, d) + min(a;, p(()) < 1. In any case, we are done. 
Now, suppose that there exists 7 < e < d such that 

{d}t < {e}+ < {h}+. (22) 

For < C < fl2i+3 a non negative integer, let e(C) = d — p((). Choose ( so that |{e}+ — {e(C)}+| 
is minimal. By (fT§)) . we know that {e(£)}+ = — {p(()} x . As moreover {d}+ — {p(a 2 i + 3)} x < 
(following Lemma 15. 2 7[) and {e(A)}+ = {h}t , we deduce that A < £ < a 2i+3 . Suppose that 
W+ - {e(C)}+ * 0. As for all C € {A,..., a2j+3 - 1}, |{e(C + l)}+ - {e(C)}+| = |{^(C)}+ - 
{fj,(C + = {q 2 i+ 2 } x , we deduce that |{e - e(C )} x \ < {q2i+2}x and the fact that \e - e(C)| < 

92i+3 contradicts the second statement of Theorem 13.251 

Thus, we have that {e}+ = {e(C)}+. Then, from we can write {e}+ = {d}+ - {fi(0}x < 
Wt = {d}t ~ {l*-W}x so that {fi(C)} x > {f-W}x- Suppose that {v(()}x > {vW} x then, as 
A < C < a2i+3, it means that C > A. But then, e = e(£) = d — /i(C) < 7 which is a contradiction 
with the hypothesis 7 < e. As a consequence, we have A = £ and e = h. 

To finish the proof, we note that (fl8|) is an immediate consequence of (fT9| and Lemma r3.27l □ 

Let a; be a rational and 7 a positive integer. From the Proposition 13 . 28l we immediately obtain 
an algorithm (see Algorithm [7} to compute the reserve ordered list of the integers q such that 
mm + (x,q)/q is a positive best approximation of x relatively to 7. 

From Algorithm [71 it is possible to obtain a bound on the number of positive best approxima- 
tions of a rational number x. In order to state the following corollary, we introduce a notation: for 
(/i, p, x) £ K 2 x N, we denote by L{pb, /?, \) the finite arithmetic sequence with first term p, common 
difference p and length \ (if X is zero then the sequence is considered as empty) . 
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Algorithm 7: Reverse order list of positive best approximations 
input 

• x = a/b = [do; a\, . . . , a n ] a rational number ; 

• the lists of integers p[k], q[k] for k = 0, . . . , n, such that are the convergents 
associated to [ao; a\ , . . . , a n ] ; 

• 7 < b a positive integer. 

output : L a reverse ordered list of the integers q such that min + (a;, q)/q is a positive best 
approximation of x relatively to 7 

1 L<- [&]; 

2 last <r- b; 

3 t <— n; 

4 if (t + 1) mod 2 = then 

5 I nextqk i — 2; 

6 else 

7 |^ nextqk i — 1; 

8 while nextqk > do 



9 

10 
11 



12 

13 
14 

15 



if last — ^[nextqk] > 7 then 

/ last — ofnextqkl — 7 
A <- floorf 



g[nextqk + 1] 
last last — A.<j[nextqk + 1] 

while last — ^[nextqk] > 7 do 
last last — q [nextqk]; 
L <- last U L ; 

nextqk <— nextqk — 2; 

16 if L[l] > 7 then 

17 |_ L ^7UL; 

is return L; 
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Corollary 3.29. Let x = [do; <2i, . . . , a n ] be a rational number, denote by Pk/qk for k = 0, . . . , n 

the associated sequence of convergents. Let 7 be a positive integer. The list a positive best approxi- 
mations of x relatively to 7 has cardinality bounded by 2 + X)!=i 2 ^ a 2i- 

Denote by L the finite sequence of increasing integers q such that min + (x,q)/q is a positive best 
approximation relatively to 7. Let L = {0, . . . , \ {ti — 1)/2J}. There exist two sequences (^j)jg/ and 
(Xi)iei with coefficients respectively in Q and N such that L = Ui G /L(/^,i, (fe+i, Xi)- Moreover, for 
i <E I , the sequence {{q}t)qeL(pi,q 2i+ i,xi) * s a ^ so an arithmetic sequence with common difference 
{q2i+i}x < 0. 

Proof. To prove the first part of the statement, it suffices to show that the number of elements 
of the list generated by the loop beginning in line [T^] of Algorithm [7] for a given value of nextqk 
is less than a next qk+i- Indeed, it is clear from the initialisation of Algorithm [7] that nextqk is 
running through the odd indices in {0, . . .,n — 1}. Now the relation g[nextqk + 1] = a nex t q k+i ■ 
g[nextqk]+(7[nextqk— 1] implies that the loop on lineQ2]is executed at most a n extqk+i times. Taking 
into account the first and last element in the list L, we obtain that its cardinality is bounded by 

The second part of the statement is clear, since the while loop on line build a (reverse 
ordered) arithmetic sequence of common difference q [nextqk] and the last point is an immediate 
consequence of (|T%|). □ 

Remark 3.30. Denote by L the output of Algorithm^ By the corollary, L is a union of arithmetic 
sequences each of which can be encoded by a triple of integers giving the first term of the sequence, 
its common difference and the number of terms of the sequence. Recall that x — [ao;ai, . . . , a n \. 
Using this encoding, the list L can be represented (as a data structure) by 0(n) bits of information. 
Moreover, it is easy to modify Algorithm [7] so that it returns the list L encoded in that way and 
have running time O(n). For this, we just have to replace lines[TW7^\ by: 



j last — 7 

length <— floor — - - 

V q|nextqk| . 

first <— last — length • q [nextqk]; 

L <— (first, q[nextqk], length) U L: 

last <— first 

We have everything at hand in order to compute efficiently the generators of JV = {x £ 
S v \v v (x) > — S/a}. Indeed, consider the line C given by the equation y + x.^ — — Let 7 — 

mod a, where 4 mod a is considered as a positive integer in {0, . . . ,a — 1}. Then —7 is the 
abscissa of the first point of the line C with integer coordinates to the left of the origin point. 
Denote by (qi)i^i the list of integers qi such that min + (/3/a, qi)j q% is a positive best approximation 
of j3 j a relatively to 7. Then if we set a, = q\ — 7, it is easily seen that the cti are precisely the 
same as the one defined in the Lemma \'S. 221 

Corollary 3.31. Let v = j3/a = [ao;a±, . . . ,a n ]. Let S be an integer. Set j¥ — {x G S u \v v (x) > 
— S/a}. Then JV is generated elements of the form {■K^ i .u ai )i e j where the cardinality of J is 
bounded by 2 + J2i2.{ 2 ^ °2i- Let I = {1, . . . , [n./2j}. There exist two sequences (/ii)i £ / and {xi)iei 
with coefficients respectively in Q and N such that (ttj)igj = Ui e jL(^, 921+1, Xi)- Moreover, the 
sequence Vv(irP i .u ai ) ai £L( fJii . q2i+1 . Xi ' ) is also an arithmetic sequence. 

By gathering all the results of this section, we obtain: 

Theorem 3.32. Let v = [ao;&i, . . . ,a n ]. Let ^# be a sub-S^ -module of S^;. Then a bound on the 
number of generators ofMax(^Sf) is rf-(2 + ^{™{ 2 ^ a,2i). These generators can be represented by d 
vectors of and d- [n/2\ arithmetic sequences of the form L(p,q,x) where q is the denominator 
of a convergent of odd index associated to [o,q; ai, . . . , a n \. 

3.3.5 Application: scalar extension of S'y-modules 

Let v' , v 6 Q such that v' > v, there is a natural inclusion 9 U ,„< : S v — > S u > . Given a module ^ over 
S y , We would like to compute the module Max(^# ®s K S u >) € Max^ ( . If M = (m,j) € M d y.k{S v ) 
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is a matrix representing it can be done by calling the algorithm MatrixReduction on the matrix 
(6 v y(mij)). 

Nevertheless, if *M is maximal, there is another better way to carry out this computation. 
Assume that ^# is represented by a couple (M' , L') with M' £ Md X k(S„) and L' = [a%, . . . , ctk] is 
a list of integers. Let (fx, . . . , fk) with /,; = w ai ■ for i = 1, . . . , k and ej £ be the basis of 
$(^) given by the column vectors associated to the couple (M',L') (see Remark |3.19[) . Then by 
definition is generated by the sub-S^-modulcs Fj = fi-S' v fl 5^. Moreover, using Algorithm [71 
one can recover a family of generators of Fj which are of the form Sj ■ ej with Sj £ and following 
Remark 13.301 it is possible to encode the generators of Fj by a list of arithmetic sequences. As this 
representation is very compact, we would like take advantage of it in order to compute the scalar 
extension. By working component by component, we only have to consider the case of a sub-SV- 
module of S v , JV — {x E S v \v v (x) > — 6/ a} for 8 £ N. Then it has been seen in Corollary 13 . 3 1 1 that 
JV is generated elements of the form (tt^ .u ai )i^j . More precisely, write v = [ao; ax, ■ ■ . , a n ] and 
let I = {1, . . . , L?i/2J}. Then, there exists three sequences (/Xj)j e j, where (Ci)ie/ an d (Xi)iei with 
coefficients respectively in Q, N and N such that (ctj)j e j — Uj e j-L(/Zj, Ci,Xi)- Let jV' — ,jV ' ®s u S v i . 
Of course, the sequence (tt^ .u aj )j £ j has coefficients in S v i and is a family of generators of J^' . 
Hence, Max(^ / ( // ) corresponds to the couple (M' , L') where the unique element of L' is given the 
minimum of all quantites /3j + v 1 ■ otj when j runs over J. Now, we remark that the sequence 
f3j + v' ■ aj is arithmeric when j runs over one subset L(fj,i, Q, Xi)- On this subset, the minimum is 
reached for the first index or the last one. Thus, to compute V , it is enough to take the minimum 
over these particular indices. It yields an algorithm whose complexity is O(n) — or 0(nd) for 
the d-dimensional case — where we recall that n is the length of the continued fraction of v (in 
particular n = 0(1 + min(log \a\, log \(3\)) if v = %,) 

3.4 Comparing the two approaches 

We have introduced two different ways to represent ^-modules and compute with them. It is 
important to compare the two approaches since they are well suited for different kind of applications. 
We call the representation of H3.2. li the (M w , M„)-representation and the representation of ij3.3l the 
(M, L)-representation. 

First, we explain how to go back and forth between the two representations. Let jtf, £ Maxg 
given with the (M, L)-presentation by the couple (M,L) with M £ Mdxk(S^) and L is a list of 
integers. We can recover a matrix Mx with coefficients in 5V whose columns vectors gives generators 
of ^# in the canonical basis of S%. Then to obtain the couple (M n ,M u ) representing we just 
have to compute the Hermite Normal Forms of Mx &>s„ S v ^ and Mx <E>s u S VtU . 

We explain how to compute the (M, ^-representation associated to a (M n) M„)-representation 
in the case that the associated module j% £ Max^ has full rank. Suppose we are given the couple 
(Mir,M u ) representing M where M v = (m nti j) £ M ( j X fc(S' i /, 7 r) and M u = (m u ,%,j) £ M dxk (S^ u ). 
We can suppose, by multiplying by a certain power of ir (which is invertible in S v ^), that all 
the m^^.j £ S v . As the coefficients of M u are defined modulo a certain power of it (namely the 
determinant of M u ), we can also suppose, by multiplying M u by a certain power of u a /n^ (which 
is invertible in S u . u ), that all the coefficients of M u belongs to S v . Let D u = det(M u ) £ S v . On 
the other side, let = det(M 3r ) / w <*.v„(d e t(M„)) g ^ By d enn ition, we have v v (D„) = 0. Denote 
by (resp. ) the sub-S'y-module of (S' v ) d generated by the column vectors of D U M V (resp. 
DttMu), considered as matrices with coefficients in S' v . We can prove: 

Lemma 3.33. Keeping the above notations, we have: 

Max((X n Jin) ® Sv S' v ) = Max(jtf£ + Jt%). 

Proof. Using the formula adj(M) = det(M).M _1 , it is clear that the column vectors of the matrix 
D U M^ (resp. D n M u ) belong to the SI ^-module generated by the column vectors of M u (resp. 
the S'^-module generated by the column vectors of M n ). As a consequence, we have JKq G 
(J( u n 'jt-n ) ® Sv S' v and Jt^ G (..#„n.4V)8s„5^. We deduce that + C (Ji u V\Ji 1 ,)®s v S' v . 
Thus, we have Max((^„ n J(^) ® Sv S' v ) D Max((^# 7r + J(£) ® Sv S' v ). 

Next, suppose that x £ Max((^f„ n J% n ) ®S„ S' v ). By Proposition 13.81 it means that there 
exists n £ N such that 7r™ • x £ (.M u n ®s u S' v and (u/vD P ) n ■ x £ (J£ u H Jt^) ®s v S' v . Note 
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that D u is a power of tt, as a consequence there exists no > n such that 

ir n ° ■ x G J% C + ..# ". (23) 

We would like to prove that there exists ni G N such that (u / w^) ni x G ^#0^ + -^o"- For this, 
it suffices to prove that (u/w^) ni x mod G Ji§ j{Jt^ R ^o") c S' v j Jt§ . As D„ is invertible 
in (remember that d„(-Dtt) = 0) there exists t € S' v and 77-2 G N such that i.ZJjr = (u/m°) n2 
mod n n °S' v . Denote by fi,...,fk the vectors whose coordinates in the canonical basis of (<S£) d 
are given by the column vectors of jMq . Now, as (u/w^) n ■ x G there exist Ai G S' uu , for 
i = 1, . . . , k, such that 

k 

But we have (u/w^) n .x G so that (u/w^) n .x G (S^) d and using the triangular form of the 
matrix M u (see Proposition I3.15[) we have that A* G 5^, for i = 1, . . . , k. By multiplying the 
preceding equation by t.D^, we obtain: 

fe 

(u/^) n+n \x + \(u/mP) n n n ° ■ x = ^(t.A,)^/,), 

i=i 

for A G S' v . Recall that we have seen that 7r"° • x G ..# 7r , thus (u/tu /3 )" + " 2 • x mod G 
./#o" / (~d?Q fl .^q). As a consequence by taking tij = n + 112, we have: 

{u/w^y 11 ■ x G jMq + -# U (24) 

By J23]) and J24|, there exists am G N such that it m -x G ..# 7r +^ " and (-u/tu' 9 ) m -a; G JZfi+JVtf. 
By applying Proposition 13.81 we deduce that x G Max(( t ^ 7r + ^#q') S^,) and we are done. □ 

Remark 3.34. In the preceding construction, we need the extension S' v of S v just to ensure that 
v v {D 7T ) — 0. Thus, if v L ,(det(M 7r )) G Z, this extension is not necessary. 

Now, let j% G Max^ be represented by a couple (M„, M u ). As and M u are given in Hermite 
Normal Form, we can easily compute and D u . Let = D U M V and M' u = D^M U . Lemma r3.33l 
tells us that we can then obtain the (M, L)-representation of jtft by calling the MatrixReduction 
algorithm on the matrix (M^M^). 

The main advantage of the (M n , M„)-representation is that is provides unique representation 
of maximal modules over S v , because of the same property for Hermite Normal Forms. Thus, it 
allows to test equality between modules. We have seen also that the echelon form is well suited to 
test whether x G Sf, is an element of ^# G Maxg as well as to computation the intersection of two 
modules. On the other side the (M, L)-representation provides an actual basis of module in Maxg . 
Moreover, the base change operation ®s v S V ' only makes sense in the (M, L)-representation and 
we will see in §Q an important application of this operation. Indeed, if v' > v, altough there is a 
natural inclusion morphism S v C S v i, the two sub-rings of § , S v>u and S u > iU are not comparable 
by the inclusion relation. 

4 Representation and precision 

In the previous sections, we have presented algorithms to compute with 5,,-modules by using, as 
a black-box, the ring operations of S v . As elements of S v can not be coded with a finite data 
structure, these procedures are not algorithms stricto sensus since they can not be implemented 
on a Turing machine for instance. In order to turn them into algorithms, we have to explain how to 
represent mathematical objects by finite data structures. Much in the same way as we compute with 
approximations of real numbers, we can represent power series with coefficients £R by truncating 
them up to a certain precision. Then we have to ensure the stability of the computations, i.e. that 
the result is independent of the part of the input that we ignore. In the following, we proceed in 
an incremental manner. First, we explain how to represent the elements of the coefficient ring $H of 
S v by a finite structure, then we deal with elements of S v and finally with more complex structures 
with coefficients in such as S^-modulcs. 
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4.1 Generality with precision 

We recall from the introduction that 9K is a complete discrete valuation ring, and that for algorith- 
mic applications we are mostly interested by: 

• Z p or more generally the ring a integer of a finite extension of Q p , 

• the ring of formal power series with coefficients in a (finite) field k. 

In any case, if tt denote the uniformizer element of d\ and p n is a positive integer, we shall represent 
an element of 9i by its image in the quotient dX/n^^K. We suppose that there exists algorithms to 
compute the arithmetic operations of the ring ^H/t^"^- We say that an element x £ D\/ir P7r 9l is 
the data of element of x £ *K up to 7r-adic precision p^ if x mod tt P7T = x. 

For the complexity analysis, we shall assume that we have efficient algorithms to perform all 
standard operations in quotients 9\/tt p *9{ for all integers p n . We discuss briefly the validity of this 
assumption for the aforementioned classical examples of rings 9\. In the case that 9\ = we 
suppose that the operations in the field k costs one unit of time and can be represented by one unit 
of memory. With that in mind, if *H = there exists a trivial algorithm to perform additions. 

It is optimal in the sense that its complexity is equal to the size of the inputs. The same thing is 
true if is the ring of integers of any finite extension of Q p . Things are more complicated for the 
multiplication of two elements of 9i/ir P7 'D\, whose time will be denoted by ToCPtt) in the rest of this 
paper. In the case *K = Z p , using Strassen algorithm [IT] , we have T{p^) — 0(p- K ) where the soft-0 
notation means that we neglect logarithmic factors. If 9\ is the ring of integer of a degree d finite 
extension of Q p , we can represent elements of 9\ with a degree d — 1 polynomial with coefficients in 
Z p and using again Strassen algorithm for polynomials, we have To(p„) = 0(d ■ p v ). If £H = 
using again Strassen algorithm for polynomials, we have T(p 7T ) = <D(p^) (we suppose here that 
operation in k costs one unit of time). We can summarize these results by saying that with the 
best known algorithms, the time To(p v ) is quasi-linear log(|lH/7r p,r lH|). 

An obvious way to obtain a finite approximation of an element of ^ aiU 1 £ is to consider 
a representative modulo a certain power p u of u. We, thus obtain a degree p u — 1 polynomial 
with coefficients in SH that we can represent by a vector of dimension p u with coefficients in up 
to precision p^ as before. We call this representation the flat approximation of an element of S v 
with u-adic precision p u and 7r-adic precision p^ or the (p u ,p 7r )-flat approximation. The data of a 
representative with 7r-adic precision p n and M-adic precision p u of an element x = ^2 a i ul /it' iul[ £ S v 
is given by a polynomial Yli=o ciiU 1 /7r^ w ~\ such that a, = a; mod n P7r . It should be remarked 
however that the flat approximation is not the only possible procedure to truncate an element of 
S u in order to obtain a finite structure. For instance, one can represent an element of S v up to 
a certain u-adic precision p u by a polynomial Y^=o 1 aiU% w i* n coefficients in SH of degree p u — 1. 
Such a polynomial may itself be represented by the data of 0.$ mod tt p ^ for i = 0, . . . ,p u — 1, as 
before but it is also possible to represent X)f=o 1 a i u% by coefficients with different 7r-adic precisions 
di mod ir p, '' i . Put in another way, we want to obtain a representative of X)f=o 1 a i ut modulo 
the *H-module X)f=o 1 7rP " i u V 71 " ' We call this representation the jagged approximation. We 
can generalize even further the flat and jagged approximations. For instance, we remark that 
for / = ^aiU 1 £ Sv the flat and jagged approximations consist in the data of f^(0)/i\ for 
i = 0, ...iPu — 1 but we could also provide the data of f^'(x)/i\ for any x £ K in the radius of 
convergence of /. 

Taking into account the previous examples, we say that a data of precision is given by any 
sub-lK-module 3? of S v . Most of the time, but not always, we want S v j 3? to be CH-module of finite 
length. Indeed, it may happen that we compute with objects of S v that can be represented exactly 
with a finite structure. This is the case for instance, if the characteristic of £R is 0, of any element 
Z C £H. In this special case, it makes sense to consider a data of precision & such that S v j is 
not of finite length in order to take into account the fact that we know certain elements of S v with 
"infinite precision". In general, in order to represent an element of by a finite data structure, 
one can consider a sub-?t- module of Sf, such that most of the time S„/ ' has finite length. 

Then, in order to compute a function / : S„ S„, we would like to replace it by its approx- 
imation / : S„l 1 — > S„l This naive approach does not work in general since, as / is not 
always £H-linear, the image by / of a data of precision is not a data of precision. Though, it is 



31 



possible to approximate f{3 2> ) by the smallest possible data of precision. One way to do this is to 
consider a regular data precision which is a data of precision that is a S^-module. Then for x £ S„ 
and one can write the first order Taylor development of / in x 

f(x + h) = f(x)+df x (h)+0(h 2 ). 

Most of the time (but not always), df x (3 2> ) will be the correct data of precision (see [T] for a full 
discussion about this). Note that any flat approximation is a regular data of precision but this is not 
always the case that a jagged approximation is a regular data of precision. The computation of the 
function / reduces to the computation of the function on the representative up to the given precision 
and the computation of the precision of the result. A more general precision data is intuitively less 
convenient for computations since it involves more complex data structures. For instance, each 
coefficient of a polynomial representing an clement of S u with the jagged approximation may have 
very unbalanced length so that it may be difficult to adapt asymptotically fast arithmetic for such 
objects. On the other side, we are going to sec shortly that even for a very common operation in 
S v such as the computation of the Euclidean division, one may take advantage of the flexibility 
of the jagged approximation. Hence, the choice of a representation to compute with elements of 
S u is a non trivial trade off between space/time complexity on the one hand and the quantity of 
precision we accept to loss on the other hand. 

It is convenient to represent a jagged precision by a series. For this, let P n = <HU % /t^- w ^ G 

S v . In the following, we denote by ^(P n ) the sub-lH-module of S v given by X^o a « u V 7r ' 
£R. Moreover, if & is sub- £K- module of S v , we denote by repr(^) : S v — > S v / ' & the canonical 
projection of £H-modules. It is clear that ^{P v ) only depends on the valuation of the coefficients 
dj of Pn = X°^q aiU l /n\- lu i € S^. If p^ is an integer, we will use the notations 3* f(p u ,Pn) for 

Pu — 1 OO 

i=0 p u 

which corresponds to the (p u ,p w )-flat approximation. If 9*' and 3* are two sub-*H-modules of 
S v such that C 3? then there is a canonical projection S v j — > S v j 3? that we denote also 
(by abuse of notation) by repr(^). If A e S u , and & is a sub- D\- module of S v , we denote by 
X.3 3 = {X.x, x 6 3 21 } the sub-$K-module of S v . If A is distinguished and S v /3? has finite length 
then S v j{\ ■ &>) has finite length. If 9>, 3*' are sub-<K-modules of S„, we denote by &> ■ 3*' the 
submodule generated by all products xy for [x, y) € {3 s x 9 s '). It is clear that if S v / 3* and S v / 3 71 ' 
have finite length then S v /{3 2> ■ 3 71 ') also have finite length. 

Lemma 4.1. For all 3 21 , 3*' sub-^K-modules of S v such that S v j and S v / s ' have finite length, 
for all x,y £ S v we have: 

1. if 3*' D 3*> then repr(^')( re P r (^')( a; )) = repr(^")(a:) / 

2. repr(^ + ^")( re P r (^)(») + repr(^ + ^")( re P r (^")(y)) = repr(«^ + + y) ; 

3. let 3* = y ■ 3* + x ■ 3*' + 3* ■ 3 s " ', then 

repr(^ > )(repr(^')(a;)) • repr(^ )(repr(^")(y)) = repr^oX^ • y); 

4. if 9" D 3*, then repr(^ 2> ')(repr(^ 2 )(a;)) • repr(^")(y) = repr(^")(a; • y) 

Proof. The fist claim is trivial. Then we have (x + 3*) + (y + 3*') = x + y + (3* + 3 s ') and 
(x + 3 s ) ■ (y + 3*') = x ■ y + x ■ 3*' + y ■ 3* + 3* ■ 3 21 ' . The fourth claim, is an immediate consequence 
of 1 and 3. □ 

We discuss briefly the complexity of the elementary arithmetic operations in S v with the (p u ,Pir)- 
flat approximation. First, we remark that the size of an element of S v with the (p tI ,p 7r )-flat 
approximation is in the order of -p u . As before, the time of an addition in S v is linear in the size 
of a representative of S v since it reduces to the addition of two polynomials of degree p u — 1 with 
coefficients in D\/ir p, 'D\. We denote by T(p u ,p 7r ) the time cost of the multiplication of two elements 
of S v with the (p M ,p T )— flat approximation. Again, by using a tweaked Strassen's algorithm, we 
have T{p u ,p 7T ) — 0(p u ■ T(p n )) — 0(p u ■ p„). In the following, we study the precision of some 
important functions using the flat and jagged approximation. 
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4.2 Finite precision computation with elements of S u 

Most of the time, even for very elementary function dealing with elements of S v , it is not possible 
to ensure the stability of the result without some extra assumption. We illustrate this fact with 
some important examples. 

First, consider the Gauss valuation function v v : K [[«]] — > Q. A natural way to define v v on 
a representative modulo & 'f(p u ,Pir), with p u ,p^ positive integers, is to compute the valuation of 
the truncated representative in S v . For instance let x — n + u 10 , then uo(repr(^ a /(9, 2))(x)) = 
vq(tt) — 1. We denote also this function by v v . But then we have uo(repr(^/(9, 2))(x)) = 1 and 
^o(repr(^ 2l /(10, 2))(x)) = 0. From the previous example, one can see that the Gauss valuation of an 
element x £ S v ^ can not be computed in general from the knowledge of its approximation. Still, it 
is possible to obtain the Gauss valuation of an element x £ S u n from the knowledge of its approxi- 
mation if we are given some extra-information about x. For instance, if v v (repr(^f(j} u ,p 7r ))(x)) = 
and if we know furthermore that x £ S v then we are sure that v u (x) = 0. More generally, it may 
happen than we have a guaranty that x £ for a A 6 Z. Then, if v is big enough, it is 

possible to compute the valuation of x from the knowledge of vepv(^f(j) Ul p 7r ))(x). 

Lemma 4.2. Let x = aiU 1 6 l/ir x ■ S v for A a positive integer. Let p u be a positive integer and 
x £ K[u] be the unique representative of x mod u Pu of degree < p u . We suppose that x ^ and 
that d = dcg w (x) < p u . 
Let is' £Q be such that 

u ,_ u >\±vM t (25) 

Pu ~ a 

then v v i (x) — v v i (x) . 

Proof. Let x = ^aiU 1 G l/ir x ■ S v . By definition, we have v v i(x) < vjc^ad) + is' ■ d and on the 
other side we have of course v v /(x) < v u '(x). Thus, in order to prove the lemma, we just have to 
check that for all i > p u , we have 

vk(cLi) + v ■ i > VK(ad) + v' ■ d. (26) 

But, using x 6 l/n x ■ S v , we get 

VK(a>i) + v ■ i > -A (27) 

for all i > p u . From (|2^|) and ([27)) . we deduce that it is enough to prove that —A + i(v' — is) > 

vxifldj + v' -d. Since is' — v > by hypothesis, it suffices to prove that —\-\-p u (y' — v) > Vk(ad) + is' .d 
for all i > p u - This is equivalent to 

, . VK{ad)+Pu-v + \ fr)Q , 

v > , (28) 

Pu - d 

which is exactly □ 

This lemma, while totally elementary, shows the following very important fact: by increasing 
the v parameter of the S^-module, one can obtain guaranties on the valuation of a certain x = 
^2 a.iU 1 £ S v from the knowledge of its representative x = X)?=i 1 a i u% with bounded Weierstrass 
degree under the general hypothesis of a lower bound on the valuation of the coefficients di. 

Another important operation for the arithmetic of S u is the inversion. 

Lemma 4.3. Let x £ S v and suppose that deg w (x) — and that v v {x) — so that by Corollary 
\2.7\ x is invertible. Let p u ,p v be positive integers. Then repr(£Pf(p Ul p 1T ))(x) £ S v / ' g f{p u ,p 7r ) is 
also invertible and we have repr(^'f(p u ,p 1T ))(x)^ 1 = repr(£?'f(p u ,p. K ))(x ). 

Proof. Write x = £ a^/Tr^ , x' 1 = £ hu 1 /^ and c = 1 = £ c^/V^ with Cj = J2l=o a i ' 
bj-i. We have vk (do) = so that we can compute ao -1 mod p v = bo mod p v . Then, using the 
formula 

bj 1 \ - , 



7rlJ v J an ^ ■nl iv \-K\ii-i>\ ' 
i=0 

together with the remark that ttI>J /(tt^J ttLGMM ) is equal to 1 or n, we obtain by induction for 
j = 1, . . . ,p u - 1, bj modp^. □ 
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Let x,y £ S u . In order to be able to compute an approximation of the Euclidean division of y 
by x , it is necessary to know that v v (y) > v v {x). One way to have that guaranty is to be given 
x with enough precision to know that it is distinguished. Then one can compute its Weierstrass 
degree. In the following proposition, we keep the notations of Proposition ^. 81 

Proposition 4.4. Let x,y £ S u . Suppose that x is distinguished and let d = deg w (x). Let q £ >S„ 
and r £ K[u] (~l S v be such that deg(r) < d and y — q ■ x + r. Put e = v v {Lo{x)) > 0, let p^ be a 
positive integer and p x = \p v je\d. Let 

Px — l OO j 

p = 7r max{p„ - \i/d\ e- \iv\ , - \iv\ }yi _l_ " 

v \iv\ ' 

p- c-d— l OO j 

p = 7r »ax{f)„ - [i/d+lj e- \iv\- \iv\ } u % , 

q ^ ^ |u/|" 

There exists an algorithm which takes as input vepr( ■^ > f(p x , p 7T )){x) and repr(^(_Pj,))(y) and oui- 
pute repr(^ , / (p x ,p T ))((7) and repr(^ , /(oo,p 7r ))(r). 




Figure 3: The form of the precision of y in the Euclidean division for d — 2 and v = 1/6. 

Proof. Recall that, from Proposition 12.81 q,r are the limits of the sequences (qj,rj) defined by 
<7o = and Tq = y and 

, Hifc.d) 
*+ 1 = *+ HiM' 



^= L °(ri)- 1 ^-Lo(»)- 



For j = 0, . . . , [Pir/el, let 



7r max{j-e- |wj |wj } yi 



P y j = 7r max{p,-Li/dJe-L"^J--L»-'J} ?i i _|_ 

i=0 i=([ P7r /el-j)-d 

and let t(j) = vepv(^(P y j))(rj). It is clear that t(0) — repr(£P(P y ))(y). We are going to prove 
that if we know repr(^f(p x ,p 7V ))(x) and t(j) then we can compute t(j + 1). Write Hi(x,d) = 
u d /it vd ■ xq, with Xo an invertible element of S„. Then from iepr(^ > f(p x ,p 7r ))(Ki(x, d)), we imme- 
diately obtain vepv(^f(p x — d,p w ))(Hi(xo, d)), and by Lemma FOI we can compute Tepi(^ l f(p x — 
d,p„))(l/(H.i(xo,d))). As &f(p x — d,p n ) C — d • j,Py,j), applying Lemma |4~T1 we deduce 

that 

vepi(^(P y<j ))(m( rj ,d)/m(xo, d)) = 

vepv(^(P yij )(Te W (^fip x ~ d, Pw ))(l/(m(x ,d)))).vepx(^(P y!j )(m( rj ,d)). (29) 
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We remark that repr(3 g (P y j)(Hi(rj, d)) = Hi(t(j), d) so that the left hand side of (j2"5)) can be com- 
puted from the known data. Dividing repr(^(P yj ))(Hi(rj, d)/Hi(xo, d)) by u^/tt^"^, we obtain 
repr (&(Pj ) ) (Hi (rj , d)/Hi(i, d) ) , where 



(r^/el-0'+l))-d-l 



7J .max{p x -L i / ci + 1 J e -L il/ J ! -L ily J u * 



i=0 



i=(rP^/el-(j'+l))'d 

Next, as u„(Lo(a;)) = e, still be applying Lemma FfTTl and remarking that t(j + 1) = ir e ■ Pj, we 
obtain 



repr(^( J + l)))(f^-Lo(x)) = 
Hi(x, a) 

repr(5»(t(j + l)))(repr(^(P J ))(Hi(r i , d)/B.i(x, rf))).repr(^ / (oo,p 7r ))(Lo(a;)). (31) 

From the above, we deduce by induction that we can compute repr(^ 2 (P a |- p7r / e ]))(r). But 
we have ^(PyJp^/e]) = ^/(oo,PTr) so that we can compute repr(^ a /(oo, p^)){r) as claimed. 
Moreover, as we can compute repr(J 2> (P ) ))(Hi(r :) , d)/Ri(x, d)), by Lemma |4?T| we can compute 
rcpr(^ ^ g (Pj))(q) — repr(Po)(g) an d we are done. □ 

Remark 4.5. In the preceding proposition, we see that in order to be able to compute repr(^'f(d,p 7r ))(r) 
we really use all the information contained in repr(&(p y , P y )){y). If we use the flat precision, then 
to obtain iej>i(^ g (d, p n )) we need to know Tepr(&([p 7r /e~\.d,p 7r ))(y). This shows that the flat preci- 
sion is not well adapted to the computation of the Euclidean division in S v since a lot of information 
about the operands is not useful for the computation. 

The proposition shows that following the computations of Algorithm Q] on representatives mod- 
ulo the given precision, we obtain the outputs with the guaranty that the result has the claimed 
precision. 

The last operation in S u (actually in S u ,ir) that we would like to consider is the gcd computation. 
To begin with, we consider some very simple examples, for elements of S v which are polynomials. 
Suppose that 9\ = Z 5 , v = so that S v = Z 5 [[«]]. Let ~F\ = repr(5 a / (oo, 2)){u - 1) and T~ 2 = 
repr(^ 2l / (cx) 1 2))(w - 2). Then it is clear that for all P 1 ,P 2 E S v such that P : = Pi mod ^/(oo, 2) 
and P2 = P2 mod ^/(oo,2) then gcd(Pi,P2) = 1. This can be seen by using the Euclidean 
algorithm to compute the extended gcd of Pi and P2 in S„/ £Pf(co, 2) which obviously returns 1. 
In this case, it is safe to claim that gcd(Pi, P2) = 1. 

Next, consider P 3 = repr(^/(oo, 2))(m — 1) and P4 = repr(^ 2 '/(oo, 2))(w — 1). In this case, 
it is very easy to find different representatives of P3 and P4 the gcd of which is not equal. For 
instance, we can take P3 = P4 = u — 1 in this case gcd(P3, P4) = u — 1 but if we take P3 = u — 1 
and P4 = u — 6 then gcd(P3,Pi) = 1. If we compute the gcd of P3 and P4 using the Euclidean 
algorithm, we obtain u — 1 and we do not have enough precision on the next remainder to decide 
whether it vanishes or not. This example shows that, in the case that the gcd of the representatives 
is not surely 1 it is not even clear how to define it since the result may change depending on the 
representatives in S v that we use in order to compute it. 



4.3 Finite precision computation with modules with coefficients in S y 

Let ^#1 and ^#2 be two maximal sub-SV-modules of Sf,. In this section, we are interested by the 
computation of the maximal sum ^#1 + max ^2 of M\ and Mi- We would like to carry out com- 
putations with finite precision and have a guaranty on the precision of the results. The preceding 
example suggests that even in the case d = 1, we can not hope much in that direction. Indeed, 
the computation of the maximal sum of two sub-S^-modules of S v ^ reduces to the computation of 
the gcd of two elements of S Uj7r and we have seen in i j4.2[ that unless this sum is S u , v , we can not 
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guaranty that the result computed with finite precision is an approximation of the result computed 
on representatives in S v ^. 

As before, we need some extra-information, that we can get from the mathematical context of 
our computation, in order to guaranty the precision of the output. A very natural extra-information 
that can arise in practise is the following: let and -#2 be two sub-S^-modules of and we 
know that there exists a positive integer c such that ^#2 C 1 /tt c ^i . We recognize a generalisation 
of the hypothesis of Lemma 14.21 where we have shown in the case that d = 1 that we can obtain 
a guaranty on the valuation v v of approximations of elements of if[[u]] for well chosen v. This 
situation is also crucial in the paper [4;. We are going to see that, although we don't know how 
to compute an approximation of ^#1 + max we can describe an algorithm which outputs an 
approximation of (.#1 ®s v S V ') +max ®s u Sv') for a well chosen v' > v. 
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Figure 4: The computation of u 1 from p u and v. 

In order to compute M\ + m ax -^2, it is enough to be able to compute + max S v ■ t where 
t £ ^#2- Indeed, let (t\, . . . , ty) be a family of generators of ^#2, we have ^\ + max ^2 = 

+max S v ■ t\ + m ax • • • +max S v • £/,,' . Let t £ ^#2 and let (e%, . . . , eh) be a family of generators of 
By our hypothesis, we know that there exists A^ £ l/n c ■ S v such that t = ^ A^. We remark 
that if all the \ are in S v then t £ jM\ so that Ji\ + S v -t — j£\ and there is nothing to do. Write 
A,; = J2j>o a ) u ^ with vk{<i)) + v ■ j > —c. Let p u a positive integer, we are going to choose i/, as 
it is explained in figure [H such that J2j> p al j u% For this it is enough to take v' > v + c/p u . 

Let t' = J^i K e i w ith A- = 2^?=o a ) u ^ ano - = Yli K e i w i tn A" = Y^p u a ) u ^ ■ Using the same 
remark as above, we have: 

{•■#1 ®S V Su> ) + m ax (t ■£„') = {J?l ® Su Sv' ) +max ■ Sy> ) +max (*" • Sy>) = {Jtl ®S V S v > ) + max (f S*), 

since t" ■ S v * £ .&\. Now, as X[ is a polynomial in u, we can obtain its valuation, greatest common 
divisor and all the operations that we need in order to compute {^M\ ®s„ Su') +max {t ■ S v >). 

We recall that we write v = f3/a with a, ft relatively prime numbers and let zu in an algebraic 
closure of K, be such that w a — ir. Let 9V = %K\m\ and = S v <S)yt 9V. The algorithm AddVector 
is an adaptation of the algorithm MatrixReduction. 

In the preceding algorithm, C'ond(X 1 L) returns true if there exists jo, Ji £ {1, ■ ■ ■ , h} such that 
Abo] • Kn] ± 0, v v (X[j ]) ^(A[ji]) - ^ and deg w (\[j }) < deg w (A[ji]) 

We want to give a consequence of this algorithm. We first need a definition. 

Definition 4.6. Lei ^# 6e a sub- S„ -module of S„. Let & be a sub-D\-module of S v . We say 
that a matrix M r = (m^) £ Mdxd'iSv/ is an £P- approximation of ^# is there exists a matrix 
M = (niij) £ Mdxd'(S u ) whose columns are the coordinates of generators of in the canonical 
basis of and such that = repr( I ^ 2 )(mjj). 

By iterating this algorithm AddVector on a set of (t\, . . . ,th>) of generators of ^#2, we obtain 
the following theorem: 

Theorem 4.7. Let ^#1 and ^#2 be two finitely generated sub-S^-modules of S 1 ^ such that ^2 C 
1/tt c ^\ for a positive integer c. Let M% — (mjj) and Mi = {mjj) be the matrices with coefficients 
in Sv of generators of ^\ and ^#2 in the canonical basis of S„. Letp Ul p^ be positive integers and 



36 



Algorithm 8: AddVector 



input 

• Me Mdxh(Su), a matrix whose column vectors C(i) for i = 1, . . . , h give generators of M\ 
in the canonical basis of S% ; 

• a list A[l], . . . , X[h] such that J2 \ l C l {M) = t, X\i] e 1/tt c • S v n and deg A[i] < p u - 1 
for i = 1, . . . , k. 

output : M e Md-Kh{S v ) and a list L a matrix such that the column vectors tu L W • Ci(M) 
give generators of + m ax t in the canonical basis of S' u 

1 L <— [0, . . . , 0]; 

2 while 3j e {1, . . . , h} such that v v {X\j]) - ^ < do 

3 while Cond(A, L) is satisfied do 

4 Pick up G {1, . . . such that A [70] ■ A [71] 7^ 0, 

^(A[io]) - ^ < ^(Abi]) ~ ^ and deg w (X\j }) < deg^Afr]); 

5 if v v (X\jo]) > v v (X\j 1 ]) then 

6 So <r- \v v (X[jo\) - w„(A[ji])]; 

7 A[j ] ^-Tr-^Abo]; 

8 L -^b'o] <~ L bo] - a • £0; 

9 (9> r ) ^— EuclideanDivision(A[7o], A[ji]); 
10 A [71] <- A [71] - qX[jo}; 

|_ C h {M)^C j0 {M) + qC h {M)- 

12 Let jo such that v v (X\j }) - ^1 = min^,...,,^]) - M) ; 

13 Let ji such that ^(A[.n]) - = min^ (A[j]) - ^zl); 

14 |_ ^b'o] <- L\j ] + av v (X\j ]) - i[jo] - a^(A[ji]) + 

15 return M, L; 
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suppose that we are given M[ = (repr(^o(P-u)P7r))(f7i^)) and M£ = (repr(£^ > o(p u ,p 7r ))(mf : j)). Let 
e be the number of columns of M2 and let v' = v + ec/p u . Then there exists a polynomial time 
algorithm in the length of the representation of M"[ and to compute a matrix = {Mf^) with 
coefficients in S v i / ^o{PmP-n) which is a &o(j> Ul P-n) -approximation of 

{Jt\ ®s„ Su>) + max {Jti ® Sv S v ,). 

Remark 4.8. If we suppose in the theorem that ^#2 is maximal, then by Theorem \ 3.3'2\ we can 
take e = d.(2 + Y)l—i a 2i) where v — [ao; ax, ■ ■ ■ , a n \. 
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